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ABSTRACT 



We expand on the results of 



Nielsen et al. 



(120081 ) ■ using the null result for 



giant extrasolar planets around the 118 target stars from the VLT NACO H and 



Ks band planet search (IMasciadri et a 



ous Differential Imaging (SDI) survey ( 



2005), the VU 



Biller et al. 



and MMT Simultane- 



Planet Survey (ILafreniere et al. 



20071 ). and the Gemini Deep 



20071 ) to set constraints on the population of 



giant extrasolar planets. Our analys is is extended to include the planet lumi- 



nosity models of 



Fortney et al. 



(120081 ). as well as the correlati o n bet ween stellar 



mass and frequency of giant planets found by iJohnson et al. 



)3) 



(120071 ). Doubling 



the sample size of FGKM stars strengthens our conclusions: a model for ex- 
tr asolar giant plan e ts wit h power-laws for mass and semi-major axis as giving 



by 



Gumming et al 



((20081) cannot, with 95% conf idence, have planets beyond 65 



AU, comp ared to the va l ue of 94 AU report ed in 



models of 



Baraffe et al 



Nielsen et al. 



Johnson et al 



mm . using the 



(120071 ) correction for 



((20031). When the 

stellar mass (which gives fewer Jupiter-mass companions to M stars with respect 
to solar-type s tars) is applie d , ho^y ever, this limit moves out to 82 AU. For the 



relatively new 



Fortney et al. 



((20081 ) models, which predict fainter planets across 
most of parameter space, these upper limits, with and without a correction for 
stellar mass, are 182 and 234 AU, respectively. 



Subject headings: stars: planetary systems 
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Introduction 



There are currentl y close to 300 extr asolar planets known, most detected by the 



radial velocity method (IButler et al 



20061 ). These planets have provided a great deal of 
information on the distribution of giant planets in short period orbits. The likelihood of 
a sta r harboring a close-in g i ant planet increase s with the metal abundance of the parent 



star (IFischer fc Valenti 



2005 



Santos et al. 



2004 ). and power laws were f ound to accurate ly 



repre sent the distributions of mass and semi-major axis of exoplanets (iCumming et al. 



20081 ). While radial velocity surveys have moved on to discovering and building up statistics 
on smaller Neptune-mass planets, direct imaging surveys continue to struggle to reach even 
the highest mass planets. 

Many observing campaigns have been conducted in the last decade to detect and 
characterize planets through direct imaging, especially aimed at young target stars, when 
the self-luminosity of hosted planets is large enough to overcome the glare of the parent star. 
Improvements in adaptive optics and instrumentation designed solely to detect planets, as 
well as specialized observing techniques, have improved the contrasts achievable close to the 
target star. This has allowed a large increase in sensitivity to planets, and resulted in the 



discovery o: 
207-3932 ( 



several planetary-rn ass (<13 Mj^p) objects, including cornpanio ns to 2MASS 



Chauvin et al. 



(jClose et al 



2007 



2004), H 



Brandeker et al. 



P 30 034 (AB Pic) 



2006 



Luhman et al 



Chauvin et al. 



20051), Oph 1622 



2007a[ ). and DH Tau flltoh et al 



20051 ). These objects were discovered with projected separations of 42, 260, 243, and 330 
AU, respectively. Few, if any, of these wide companions (>200 AU for objects found around 
stars) are likely to have formed in the primordial circumstellar disk of the primary. 



Recently there have been exciting discoveries of planetary-mass ob 



2008 



ects around the 



Kalas et al. 



2008 



higher-mass A stars: HR 8799, Fomalhaut, and B Pic ( Marois et al. 
Lagrange et al.l 120091 ). In the case of the triple planet system HR 8799, we know that the 
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largest separation planet 



parallax as the primary (jClose fc Males 



(HR 8799 b, with a 



projected separation of 68 AU) has the same 



20091) and so it (and very likely HR 8799 c and d) 



formed together around the A5 star HR 8799. While this is an amazing system, in this 
study we concentrate on lower mass stars, more similar to the Sun. 



In 



Nielsen et al. 



([20081), we presented nuL 



for extrasolar giant planets of 



Masciadri et al 



results fro m the direct imagi ng surveys 



fl2005h and 



Biller et al. 



(120071 ). using the 



contrast curves for each of 60 unique target stars to set constraints on the populations of 
extrasolar planets. We concluded that extrasolar gia nt planets are rar e at la rge separations 



Lafreniere et al. 



(120071 ) published the 



(>60 AU). Just prior to the publication of our work, 
null results from the Gemini Deep Planet Survey (GDPS) for 85 stars, reaching conclusions 
very similar to ours. In this paper, we combine the samples from these three surveys, to 
improve the statistical constraints we can place on extrasolar giant planet populations. 

Also since these past publications, two i mportant pape r s hav e been published relevant 



to direct imaging of extrasolar giant planets. 



Johnson et al. 



(120071 ) compared radial velocity 



target stars of different masses, and found that less massive stars have a lower likelihood of 
hosting a giant planet (>0.8 Mj„p). Since direct imaging surveys lean heavily on the M 
stars in their samples (as these are intrinsically fainter, making the detection of close-in 
planets easier), we attempt here to estimate the corresponding decrease in the strength of 
earlier null results. Also, a new set of planet luminosity models have bee n published by 



Fortnev et al 



(120031) and 



(120081). wh i ch dif fer from the popular "hot start" models of 



Baraffe et al. 



Burrows et al 



(120031 ) which had been previously utilized in such work. These 



new models are based heavily on the "core accretion" model of planet formation, and tend 
to predict consistently fainter fluxes for giant planets, especially at the youngest ages and 
largest planet masses. In addition to significantly enlarging the sample, this paper takes 
into account the stellar mass dependence of planet frequency, and the new core accretion 



models of iFortney et al.l (120081 ). to present more realistic constraints on the distribution of 
extrasolar giant planets around Sun-like stars. 



Observations 



2.1. VLT NACO H and Ks imaging 



Masciadri et al. 



(120051 ) carried out a survey of 28 young, nearby, late-type stars with 
the NACO adaptive optics system at the 8.2 meter Very Large Telescope (VLT). These 
observations have exposure times of order 30 minutes, with stars being observed in the 



H or Ks bands. 



Masciadri et al. 



or t he 22 stars used (see Section 12. 4p from the VLT NACO survey of 



(120051 ). the median target star is a 12 Myr old K7 star at 30 pc. 



2.2. VLT NACO and MMT SDI 



A survey of 54 young, nearby stars of a variety of spectral type s ('between A and 



M) 



was conducted between 2003 and 2005, with the results reported in iBiller et al.l (120071 ). 
This second survey used the Simultaneous Differential Imager (SDI) at the 6.5 meter 
MMT and the 8 meter VLT, an adaptive optics observational mode that allows higher 
contrasts by imaging simultaneously in narrow wavelength regions surroundin g the 1.6 Atm 



methane feature seen in c ool brown dwarfs an d expected in extrasolar planets (ILenzen et al. 



2004 



Close et al. 



20051 ) (ISwain et al 



(120081 ) have recently detected methane in the 
atmosphere of the transiting extrasolar planet HD 189733b). This allows the light from 
a hypothetical companion planet to be more easily distinguishable from the speckle noise 
floor (uncorrected starlight), as the two will have very different spectral signatures in this 
region. This translate s to higher sensitivity at smaller separations than the observations of 



Masciadri et al. 



(120051 ). which were conducted before the VLT SDI device was commissioned 
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(see Fig. 14 of iBiller et al.l (l2007l ) for a more detailed comparison of the two surveys). For 
most of these SDI targets, the star was observed for a total of 40 minutes of integration 
time, which includes a 33 degree roll in the telescope's rotation angle, in order to separate 



super speckles-which are created within the instr ument, and so wi 



physical companion, which will rotate on the sky (iBiller et al 



1 not rotate-from a 



from the 



Biller et al. 



20061). The 50 stars used 



(120071 ) SDI survey have a median age, distance, and spectral type of 



70 Myr, 24 pc, and Kl, respectively. 



2.3. Gemini Deep Planet Survey 

At about the same time as the 



Biller et al. 



(120071 ) SDI survey, a direct imaging 



campaign was underway from the Gemini North telescope using the Altair AO system 
and NIRI camera, imaging in a narrow-band H filter with transmission between 1.54-1.65 
Hm. The observations were done using the Angular Differential Imaging (ADI) technique, 
which leaves the Cassegrain instrument rotator off during a sequence of exposures on 
the star, so that instrumental effects like super speckles will stay fixed, w^ 



companions (like planets) w ill rotate throughout the observation (iLiu 



2006 



Lafreniere et al. 



fixed, wh 


ile physical 


(Liu 


2004; 


Marois et al. 



20071 ). This technique is most effective at producing high contrasts 



as one moves awa y from the star, with the contrasts achieved exceeding those with SDI 



Biller et al 



Lafreniere et al. 



2007) bey ond ~0.7". For the 71 of the 85 stars from the GDPS survey of 



(120071 ) which we consider here, the median target star is a KO star at 
a distance of 22 pc, with an age of 248 Myr. Hence, the target stars of this survey are 
somewhat closer and older, whereas the southern VLT SDI survey was focused on more 
distant, though younger, stars. 
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2.4. Target Stars 



Between the three surveys hsted above, our analysis considers 118 distinct target stars, 
with some overlap between surveys. General properties of the target stars, including name, 
position, distance, spectral type, age, fluxes, and observation method are given in Table [TJ 
We a ttempt to derive age s in a uniform manner for all target stars, using the same method 



as m 



Nielsen et al. 



(120081 ). If the star is a member of a known moving group, the age of 
that group is adopted as the age of the star. We limit moving group identiflcations to 
the well-studied and established groups AB Dor, Her/Lyr, Tuc/Hor, jS Pic, and TW Hya. 
Mem bership in more c ontro versial associations, such as the Local Association and IC 2391 



[e.g. 



Fernandez et al. 



(120081 )). are not adopted here. If the star is not a member of a group, 
but has a measured value of the calcium emission indicator K^hk and a measurement of the 
equivalent width of the lithium absorption at 6708 A, the average of the ages from the two 
methods is used. If only one of these two spectral age indicators is available, the age from 
that measurement is used. If a star from any of the three surveys has none of these three 
sources for an a g e esti mate, it is simply not us ed in this w o rk. A s a result, 6 stars from the 



Masciadri et al. 



the 



Lafreniere et al. 



(|2005r) sur vev. 1 star from the 



Biller et al. 



(120071 ) survey, and 14 stars from 



(120071 ) survey were dropped. 



by 



In order to determine ages fro m the H^hk value, we utilize the polynomial flt derived 



Mamajek &: Hillenbrandl ( l2008l ). The authors derive their relation from K^hk values of 



young clusters, and flnd a precision of 0.2 dex for ages derived from this relation. 

For lithium values, we compare the equivalent width of the 6707 A lithium line and 
effective temp erature of the star to a set of young stellar clusters. For each clu ster (NGC 



2264 - 3 M vr (ISoderblom et al 



125 Mv r (ISoderblom et a. 



300 Myr (ISoderblom et al. 



19991), IC 2602 - 50 Mvr ( 



1993ah . M34 - 250 Mvr (I Jones et a 



Randich et al. 



1993bh . M67 - 5200 Myr (l.Tones et al 



200 ih . Pleiades 



19971) . Ursa Majoris 



19991 )). the mean lithium 
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equivalent width is fit as a function of effective temperature. Tlien, for our target stars, we 
interpolate between tlie fits to eacli cluster for that star's effective temperature, and the 
lithium value gives us the age (E. Mamajek private communication). 

For target stars with both a lithium and an K^hk age measurement, the median scatter 
between the two is a factor of 3. When we consider stars in our target list that belong 
to a single moving group (e.g. AB Dor or /3 Pic), and compute their ages using only the 
lithium or K^hk method (that is, we temporarily ignore their membership in a group), we 
find the scatter in the computed age, between members of the same moving group, to also 
be about a factor of 3. This suggests that the noise in our age measurements is primarily 
astrophysical in nature. While finding a precise age for any single target star is notoriously 
difficult, our hope is that by using a large sample of stars the individual errors will average 
out of our final results. 

Table [2] gives details on measurements (if available) for each of the three age 
determination methods used here, as well as the final adopted age, for each target star. We 
also plot our targets in Fig. [T], giving the age, distance, and spectral type (using absolute 
H magnitude as a proxy) for each star. Overall, for all 118 of the stars considered in this 
paper, the median target star is a Kl star at a distance of 24 pc with an age of 160 Myr. 



3. Monte Carlo Simulations 



As in 



Nielsen et al. 



(120081 ) we use Monte Carlo simulations of "fake" planets around 
each of the target stars in the three direct imaging surveys considered here. A large number 
(10^ - 10^, depending on the application) of simulated planets are given random values 
of eccentricity, viewing angles, and orbital phase based on the appropriate distributions. 
Planet mass and semi-major axis are assigned either from a grid (see Section [3.31) . or 
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from power- law distributions (as in Section [33]). For graphical representat ions of the 



Nielsen et al. 



distrib utions of extrasolar planet orbital parameters, see Fig. 2, 5, and 6 of 
( 120081 ). and the discussion therein. For each observation of a given target star, the flux 
of each simulated planet is computed based on the planet's mass and the target star's 
age, using one of three planet models (see Section ISTTIl . The angular separation between 
parent star and simulated planet, as well as the flux ratio between planet and star, are then 
computed given the distance to the star. These are compared to the contrast curve for the 
observation, which give the faintest detectable companion (at the 5a level) to the star, in 
the observation band, as a function of angular separation from the star. 

In cases where the same target star is observed in multiple epochs, and sometimes 
among different surveys (a common occurrence, the 22, 50, and 71 stars we use in this 
work would suggest a sample size of 143 target stars, but there are only 118 unique target 
stars between these three surveys with reliable age estimates), the additional elapsed time 
is taken into account. Simulated planets are generated at the earliest epoch as usual, and 
compared to that contrast curve. Their parameters are then used again, with orbital phase 
advanced forward by the time between observations (often a small effect for the planets to 
which these surveys are sensitive, a 30 AU orbit around a solar-type star has a 160 year 
period, and the typical time span between observations is at most about 3 years), the fluxes 
of the simulated planets are now computed in the new observation band, and compared to 
the new contrast curve. The process is repeated for as many observations as were conducted 



of the target star, and a simulated pla net that is detectab 



epochs is considered detectable. Again, 



Nielsen et al. 



e in any of the observational 



(120081 ) provides more details on these 



simulations, in particular their Fig. 3, 4, and 7. 
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3.1. Theoretical Models of Giant Planet Fluxes 



In order to use the measured contrast curves for each observed target stars to determine 
which simulated planets could be detected, it is necess ary to have a conversion from planet 
mass and age t o NIR flux . As in 



Burrows et al. 



(20031) and 



Nielsen et al 



Baraffe et al. 



(120081 ) . we use the theoretical models of 



(120031 ) for the calculation of exoplanet flux, using 



the mass of each simulated planet and the age of the host target star, using the fllter 
band (H or Ks) appro priate for t he particular obs ervation. In the cases of the GDPS 



(ILafreniere et al 



20071 ) and SDI flBiller et al 



20071 ) surveys, where the observation band 



was a specialized filter instead of the standard H bandpass, a correction factor is applied 
(see Section [3.21 for details). Though these two "hot start" models provide basically 
similar predictions, we perform our calculations with both, as the two models can predict 
significantly different NIR fluxes for exoplanets, depending on planet mass and stellar age, 
as shown in Fig. [2J 



Since the publicati on of 



have been published by 



Nielsen et al 



Fortney et al 



j2008l ) . an additional set of theoretical models 



(120081 ) for extrasolar planets for a r ange of masses 



and a ges 



fl2003h and 



he maior difference between these new mode^ 



Baraffe et al. 



(120031 ) is that the 



Fortney et al 



s and those from iBurrows et al. 



(120081) mode l s are based heavily 



Hubickyj et al. 



(120051 ) )■ where giant 



on the "core accretion" theory of planet formation (e.g. 
planets are formed from an initial ~10 core accreting gas from the protoplanetary disk. 
After the brief luminous accretion phase, these models predict consistently fainter NIR 



fluxes than the "hot start" models 
nicely, see Fig. 1 of iFortney et al. 



unti 



■'lOO Myr to ~1 Gyr, when the models overlap 



( 120081 )). which do not base their initial conditions on 



planetary core accre tion models. For more detail, consult Figure 8, and Tables 1 and 2, of 



Fortney et al. 



(120081 ). 



I 

As is the case with the iBurrows et al 



(I2OO3I ) models, the iFortney et al. 



(120081) models 
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do not cover the full range of planet param eters we consider her e (masses between 0.5 and 
15 Mjup, ages from 1 Myr to 10 Gyr), since 



Fortney et al. 



(120081 ) limit their calculations to 



planets with Teff >400K, leaving the consideration of cooler planets to future work. As a 
result, we extrapolate the models to masses below 1 Mj„p and above 10 Mjup, and at larger 
ages (the age a planet cools below 400 K depends on the mass of the planet, ~30 Myr for 
a 1 Mjup planet, and ~1 Gyr for a 10 Mj„p planet). While not an ideal solution, as we 
are ignoring the complicated physical processes taking places in planets as we cross these 
boundaries in exchange for simple relationships between NIR fluxes and age and mass, we 



believe that overall this met 
planets as predicted by the 



lod provides a good overall picture of the fluxes of extrasolar 



gridpoints of both the 



Fortney et al. 



Baraffe et al 



(200 



(120031 ) and 



models. In Fig. [2l we plot the initial 



Burrows et al. 



( I2OO3I ) models, as well as our 



extrapolation s to the full range of par ameter space. A simi lar plot comparing the predicted 



fluxes for the iBaraffe et al. 



(l2003f ) and 



Fortney et al 



(I2OO8I ) models is shown in Fig. [31 Our 



effort to map additional areas of model parameter space is worthwhile since this work is the 
first to apply these new core accretion models to the field of high contrast imaging surveys. 



3.2. Narrowband to Broadband Colors 



Nielsen et al 



(120081), we 



When we considered stars observed with the SDI method in 
used a constant conversion from the broadband H magnitude predicted by the models 
to the measured contr ast in the narrow band "off-methane" filter (SDI Fl, 2% bandpass. 



centered at 1.575 /im 



observed T6 objects (IBiller et al 



Close et a. 



20051) ). While this conversion factor was consistent with 



20071 ). it would be expected to vary across a broad range 



of planet temperatures, corresponding to the large differences in ages and masses of the 
simulated planets. In this work, we used template spectra from the SpeX instrument of 
132 low-mass objects, spanning spectral types from LO to T8, to compute the difference 
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between broadband H and narrowband filters as a function of effective temperature (M. 
Liu, private com munication). Spec t ral ty pes are converted to effective temperature by the 



polynomial fit of 



Golimowski et al. 



(120041 ). the ir Table 4. S peX s p ectra were obtained from 



the online SpeX Pris m Spectral Libraries (e.g. 



Cruz et al. 



(120041 ) 



Kirkpatrick et al. 



(120061 ) 



and iBurgasserl ( l2007l )). Since the reliability of the models at reproducing the methane band 
when modeling planet atmospheres is still uncertain, we prefer this method to purely using 
the synthetic spectra from the models to make this color correction. 

For GDPS target stars, we use this conversion for the NIRl CH4-short filter, to convert 
the model's prediction of planetary H-band flux to this 6.5% bandpass filter, centered at 
1.58 /im. For SDl target stars, we follow the steps of the data reduction used in computing 
contrast curves, computing planet fluxes for both the bluest "off-methane filter" (Fl), and 
the "on methane filter" (F3) both with a 2% bandpass, centered at 1.575 and 1.625 fim, 
respectively. Just as is done for the survey images, the on-methane flux is subtracted from 
the off-methane flux, providing (for each value of effective temperature ) the expect e d fina l 



Biller et al. 



(120071 ). 



flux in the subtracted image, as represented by the contrasts curves of 
For both the SDI and GDPS target stars, we use the appropriate effective temperatures 
predicted by the models to match these color corrections to simulated planets of each 
combination of age and mass. 



To partially account for this effect in 



Nielsen et al. 



(120081 ) ■ for SDI target stars, we 



had imposed an upper cut-off on planet mass, set by where the models predicted planet 
effective temperatures would rise above 1400 K for a given age. Above this temperature, 
the methane break would be so weak that subtracting the "on-methane" image from the 
"off- methane" image would simply remove all flux from the planet, as it is meant to do 
for the star. As a result, planets more massive than this limit were simply considered 
undetectable. With our more robust method that appropriately attenuates planet flux as a 
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function of temperature, it is no longer necessary for us to impose this rather crude binary 
cut for SDI targets. 

In principle, an SDI observation of a non-methanated companion should not suffer 
from self-subtraction of the companion signal, as images in the three SDI filters are shifted 
by wavelength before subtraction. This step aligns the speckles in the images (which scale 
as ^, where A is the observation wavelength and D is the diameter of the telescope), but 
misaligns any physical objects (where separation from the primary star on the detector is 
not a function of wavelength). As such, following subtraction of images from two different 
filters, a real companion should appear as a "dipole:" a positive and negative PSF, forming 
a radial line toward the primary star. The separation between the positive and negative 
parts of the dipole in the subtracted image would be given by Ad ~ where Ad is the 
length of the dipole on the detector, AA is the difference in wavelength between the two 
filters, and d is the separation on the detector between the primary star and the companion. 
The most extreme shift in filters for SDI observations is between the 1.575 fim and 1.625 
fim filters, or 3%. Since the field of view for the NACO VLT SDI observations was only 
2.5", the largest shift between positive and negative companions in the subtracted image 
would be 6.5 pixels. As the FWHM for these observations was typically 3.5 pixels, this 
dipole effect would easily be lost against the speckle background at large separations, and 
almost undetectable at small separatio ns, where positive and negative companions would 



more closely overlap (IBiller et al 



20071 ) 



3.3. Completeness Plots 



Using a similar method to iNielsen et al.l (l2008l ). we run Monte Carlo simulations of 
extrasolar planets at a grid of mass and semi-major axis points for each target star. For 
each star, then, we have what fraction of simulated planets could be detected as a function 
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of planet mass and semi-major axis. In order to combine these results over all 118 target 
stars, we again make use of the concept of the "planet fraction," or fraction of stars with a 
particular type of planet, defined such that 

N{a,M)= Yl fpia,M)P,{a,M) (1) 

i=l 

where N(a,M) is the number of planets we would expect to detect, as a function of 
semi-major axis and planet mass. Nobs is the number of stars observed, and Pj(a, M) is the 
fraction of simulated planets, at a given combination of planet mass and semi-major axis, 
we could detected around the ith star in the sample. fp{a, M), then, is the fraction of stars 
that have a planet with a mass M and semi-major axis a. If every star had one Jupiter-mass 
planet at 5 AU, for example, then fp{5AU, IMjup) = 1, and the number of these Jupiter 
analogs we would expect to detect from the three surveys would simply be the sum of the 
detection efficiency for these planets around all target stars. That is, if we had 10 stars in 
our sample {Nobs = 10), and we had a 50% chance of detecting a Jupiter-like planet around 
each star {Pi{5AU, IMj^p) = 0.5), our expected number of detections of these planets would 
be 5. 

In the case of not finding planets, as was the case for the three surveys of FGKM stars 
considered here, we can use the null result to set an upper limit on the planet fraction, 
fp. If we assume that planet fraction is constant across all stars in our survey (we will 
reexamine this assumption in Section 13.51) , we can remove fp from the sum of Equation [H 
Then, utilizing the Poisson distribution, where the probability of detections given an 
expectation value of 3 (that is, N{a,M) = 3), is 5%, we can set the 95% confidence level 
upper limit on planet fraction with the equation 
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fp{a,M) < 



(2) 



So, with the above example, where the expectation value is 5 for Jupiter-analogs over 10 
target stars {J2i'=i~^^ Pi{^^U,lMjup) = 5), not detecting any such planets would allow 
us to place a 95% confidence level upper limit of 60% on the fraction of stars with a 
Jupiter-twin (/p(a,M) < |). Doing this over the entire grid of planet mass and semi-major 
axis allows us to plot what constraints can be placed on combinations of these planet 
parameters. 

Fig. m gives the upper limit on pl anet fraction a s a fu nction of planet mass and orbital 



semi-major axis, usin g the models o f 
theoretical models of Burrows et al. 



Baraffe et al 



(120031 ) . with a similar plot using the 



(120031 ) given in Fig. O We can place our strongest 
constraints on planets more massive than ~4 Mj„p between 20 and 300 AU (fewer than 
5% of stars can have such planets at 68% confidence); when stars of all spectral types 
are considered, the lower limit probed by direct imaging and the upper limit of the radial 
velocity method are still a fac tor of 5 apart. When we repeat the calculations using the 



models of 



Fortney et al. 



(120081 ) ■ the decreased NIR fiux predicted for giant planets reduces 
constraints that can be placed on extrasolar planets, with the "sweet spot" moving out to 
~80 AU, as seen in Fig. O 

While we continue to run calculations using all three sets of models, and report the 



results here, 



Baraffe et al. 



or the sake of brevity we will henceforth only plot figures corres 



(12003[ ) C QND models. 



jjonding 



to the 



mom "hot-start" and 



Jowever, the figures appropriate to the iBurrows et al. 



Fortney et al. 



(120081 ) core accretion models are available in our 



supplement, available at this URL: http:/ /exoplanet. as. arizona.edu/~lclose/exoplanet2. html 
The supplement also contains individual completeness plots for each of our 118 target stars, 
using each of the three models of planet fiuxes. Additionally, we summarize basic results 
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for all of our calculations in Table [31 



3.4. Testing Power Law Distributions for Extrasolar Planet Mass and 

Semi-Major Axis 

These null results for extrasolar planets are also useful in setting constraints on the 
parameters of models for planet populati ons that assume power law distributions for the 
semi-major axis and mass distributions. 



Gumming et al. 



(120081 ) carefully examined the 



sensitivity of the Keck Planet Search, and determined that, over the range to which the 
radial velocity technique is sensitive (0.3 to 10 Mj^p, 2-2000 day orbital periods), planets 
follow a double power-law distribution with index -1.31 in mass and -0.61 in semi-major axis 
(-0.74 in orbital period). That is, oc M~^'^^ and 4^ oc a~°'^^ (note t hat we define pow er 



l aw in dices with re spect to linear b i ns, not the logarithmic bins of iCumming et al 



(120081). Also, while 



Gumming et al. 



(120081 ) use a and (3 to refer to the power law indices for 



mass and period, respectively, we use a to refer to the power law index for semi- major axis). 



Binarity is likely to disrupt planet formation, or at the very least 



change the underlying 



Bonavita &: Desidera 



distrib ution of planets between binary planet hosts and single stars. 

(120071 ) have shown that the distribution of radia l veloc ity planets for binary and single-star 



hosts are quite similar, and 



Holman fc Wiegert 



(I1999I ) suggest that planets are stable in 



binary systems with a planet semi-major axis >20% of the binary separation. We take 
this into account for our consideration of power-law distributions of semi- major axis by 
excluding target stars with binaries within a factor of 5 of the planetary semi-major axis 
being considered. In Table HI we give the results of a literature search for binaries among 
our target stars, including binary separation and binary type. 



By adopting these power laws, and using the normalization of 



Fischer fc Valentil (120051 ) 
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to give the total fraction of stars with planets, we can then predict how many planets these 
three surveys should have detected for various power law fits. If a large number of planets 
is predicted, our null result can be us ed to strongly exclude that model. If we accept 



the power-law dis t ributi on for mass of 



Gumming et al. 



(120081 ) and the normalization of 



Fischer fc Valentil (120051), the two remaining parameters are the semi-major axis power-law 
index a, and the semi-major axis upper cut-off (that is, what maximum semi-major axis the 
distribution continues to until planets are no longer present). We illustrate this in Fig. [3, 
where we depict various models of the semi-major axis distribution, an d the confidence 
with which we can reject them, using the models of 



Baraffe et al. 



(120031). For 12 different 



combinations of semi-major axis power law index and upper cut-off we give the percenta ge 



we ca n reject each of these 12 models in this figure. For the model of iCumming et al. 
( 120081 ) ■ with ^ oc a~°-^^, and at 95% confidence, the upper cut-off must be less than 65 
AU, and less than 30 AU with 68% confidence. 



We again use the theoretical models for planet fluxes of 



Baraffe et al. 



(120031), and 



consider a broader range of power-la w index a and upper cut-off in Fig. [HI A s before, the 



2003 



Baraffe et al. 



20031 ) are generally 



results from the two hot start models (iBurrows et al. 
similar, as the upper cut-off s must be less than 28 and 56 AU (68% and 95 % confidence) 



for the 



Burrows et al 



( I2OO3I ) models. The fainter predicted fluxes from the 



Fortney et al. 



( 120081 ) models reducing the areas of parameter space that our null result can exclude: the 
68% and 95% confidence level upper limits for upper cut-off become 83 and 182 AU, for a 
-0.61 power law index. 
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3.5. The Dependence on Stellar Mass of the Frequency of Extrasolar Giant 

Planets 



In Sections 13.31 and 13.41 we assume t he distribution and fr equency of giant planets is 
constant across all the stars in our survey. 



Johnson et al 



(120071 ) show this assumption to be 



incorrect by examining the frequency of giant planets around stars in three mass bins from 
radial velocity surveys, and show ing that more massive stars are more likely to harbor giant 



planets (see their Fig. 6). As in iNielsen et al.l (120081 ). we divide the target stars into two 



samples, one containing only M stars, and the other with FGK stars (Our sample contains 
a single A star, HD 172555 A, with spectral type A5, with all our stars F2 or earlier. We 
include this A star with the FGK stars; however, observations of this single star are not 
sufficient to make any meaningful statements about the population of planets around A 
stars). We then imagine a planet fraction (/p) with one value for M stars, and another for 
stars of earlier spectral types. 



In Fig. [To] we use the iBaraffe et al.l (120031 ) models to show the upper limit that can be 
placed on the planet fraction for M stars. Since only 18 of the 118 target stars are M stars, 
the smaller sample size greatly reduces the constraints that can be placed on planet fraction 
near the center of the contours, and the outer edge in semi-major axis. Interestingly, 
the small separation edge of the contours is virtually unchanged between Figs. H] and 
[TOl indicating that the power with which these surveys can speak to the populations of 
short-period giant planets is entirely due to the M stars in the surveys. 



Fig. [TT] uses the models of 



Baraffe et al. 



(I2OO3I ) to give the upper limit on planet 



fraction for the FGK stars in the survey. The result for long-period planets and within the 
central contours is much the same as for stars of all spectral types (Fig. H]), but the contours 
at the smallest values of semi-major axis march outward without the M stars to provide 
high contrasts at small angular separations. 
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A more satisfying way to address the issue of stellar mass dependence is to weight 
the results by target star mass, so that all stars in the survey can be applied to the 
result simultaneously. To do this, we construc t a linear fit to the metallicity-corrected 



histogram from Fig. 6 of 
a function of stel 



Johnson et al. 



by 



Johnson et al. 



(j2007l ). to give a correction to planet fraction as 
ar ma ss, as we show in Fig. [9l (Here we assume that the relation found 
(120071 ) for short-period planets (less than six years) applies to the entire 
range of semi-major axis. While this assumption is obviously untested, in the absence of 
better data we believe it is a good starting point.) In the case of setting upper limits on 
planet fraction, we now allow planet fraction to become a function of stellar mass (M*) in 
addition to planet mass and semi-major axis {Mp and a). In that case we can specify planet 
fraction for the stellar mass of a solar mass (/p,i.o(fl, Mp)), and find the upper limit as with 
Equation [2], but now including an extra term for the mass correction: 



/„ 1 o(a, Mp) < — ^ (3) 

' " Ef=i P^i(^, Mp)mcMM.,) 

where mci.o(M* j) is the mass correction as a function of the stellar mass of the ith star 
in the sum, normalized to 1.0 Mq, and defined by mci.o(M*) = -p^^^^^, where Fp is the 
fraction of stars wi th a detected r a dial v elocity planet as a function of stellar mass, using 



the linear fit to the 



Johnson et al. 



(120071 ) results. Again, going back to our earlier example. 



imagine that we have 10 stars, each with 50% completeness to Jupiter-like planets. If all 
10 stars are 1 solar mass, then mci,o{M^) = = 1, and as before the upper limit 

on planet fraction (for the 95% confidence level, as given by the 3 in the numerator) is 
60%. On the other hand, if only four of the ten target stars had masses of 1 Mq, and the 
remaining six had masses of 2.5 Mq, we must weight the results to account for the greater 
likelihood of stars of earlier spectral types to have planets. A 2.5 Mq star is twice as likely 
to have a planet as a solar mass star (see Fig. [9]), so for the four stars of 1 Mq, mci,o{1.0MQ) 
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remains 1, as before, while for the stars of 2.5 Mq, this factor doubles, mci.o(2.5M0) = 2. 
In this case. Equation [2] becomes fp,i.o{a,Mp) < 0.5+0.5+0.5+0.5+1+1+1+1+1+1 = 3/8 =37.5%. 
Including A stars in this fictional example almost doubles the constraint we can place on 
the fraction of stars with a giant planet. Similarly, M stars will be weighted against to 
account for their decreased likelihood of having planets. As an aside, we note that while 
our sample is spread across spectral type (1 A star, 8 F, 33 G, 58 K, and 18 M stars), only 
18 of our target stars are more massive than the sun. Despite the increased probability 
of finding planets around higher mass target stars, these stars are intrinsically brighter, 
and so moving earlier in spectral type very quickly results in any potential planet photons 
being swamp ed by the glare of it s host star (though the recent discoveries of planets around 



A stars, e.g. 



Marois et al. 



(120081 ). show that this difficulty can be overcome and produce 



exciting results). 

In Fig. [T2|, we plot the upper limit on planet fraction for stars of 1 Mq, using 
Equation [3l When comparing Figs. [10] and [11] with Fig. |U we see that the contours at small 
values of semi-major axis are set mainly by the 18 M stars in our sample, while the behavior 
at large separations and the depth of the contours at intermediate values of semi-major axis 
are set by the 100 FGK stars in the sample. So it is then not too surprising that Fig. [12] 
is quite similar to Fig. H] with the contours corresponding to the smallest upper limits on 
planet fraction shrinking slightly, and the contour at lower semi-major axis moving to the 
right in the figure, as M stars are now given less weight. 

Alternatively, instead of normalizing to solar-type stars, we can instead consider what 
constraints are placed on stars of 0.5 M© (about an MO spectral type). The constraints 
should become more powerful, as we assume a global decrease in the planet fraction for 
massive planets around lower-mass stars. (Again, this applies strictly to massive planets, 
>0.5Mjup. The direct imaging surveys considered here are not sensitive to Neptune mass 
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planets, which may be more common around M-stars: lEndl et al.l (120081 ) suggest that Hot 
Neptunes may be ~4 times more prevalent orbiting M-stars than Hot Jupiters around 
FGK stars) In fact, the only result of this change is to multiply a constant factor by the 
right-hand-side of Equation [3] corresponding to the ratio of the likelihood of finding a planet 
around a solar mass star to that of finding a planet around a star of 0.5 M0, or 1.5 in this 
case. We plot these lim i ts on planet fraction for half solar mass stars in Fig. [T31 with the 



models of 



Baraffe et al 



(I2OO3I ). As expected, the contours move outward, setting strong 



constraints on the frequency of giant planets in long-period orbits around M stars. 

We also reconsider the implications of stellar mass on the constraints put on the 
power-law model for the semi-major axis distribution of e xtrasolar planet s , as d iscussed in 



Section 13.41 Again, by using the linear fit to the results of 



Johnson et al 



(120071 ). we boost 



the predicted number of planets for higher mass target stars, and suppress that number 



three power 



Baraffe et al. 



aw indices 



20031), but 



for lower mass stars. In Fig. [H] we show the same combination of 
and four values of the upper cut-off as before, and the models of 
now with the additional correction for the dependence of planet frequency on the stellar 
mass of each target star. The confidence level at which we can exclude each model drops 
compared to Fig. [71 as M stars are effectively given less weight. While the upper limit on 
planet fraction as a function of planet mass and semi-major axis is specific to a given stellar 



mas s. Fig. (and t 



The 



Johnson et al. 



le ne xt one. Fig. [T5|) need not 



3e normalized to a specifi c spectral type. 



Fischer fc Valentil (120051 ) planet fraction 



(120071 ) mass correction and the 
sets the absolute likelihood a given target star has a giant planet, which is used to calculate 
the predicted number of planets detected from our entire survey. Given our null result, 
this expectation value is used to set a confidence level with which the entire model (giant 
planet self-luminosity, giant planet fraction, dependence of planet fraction on stellar mass, 
and planet mass, semi-major axis, and orbital eccentricity distributions) can be rejected. 
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For the full range of power law index and upper cut-off, again with the 



Baraffe et al. 



( 120031 ) models, we plot contours for the confidence level of rejection in Fig. [T3 This figure 
is again generally similar to Fig. [HI b ut with the con s traint s slightly looser as M stars in the 
sample receive less weight. With the 



Johnson et al. 



(120071 ) mass correction, the 68% and 



95% confid ence level upp e r limi ts on the semi-major axis distribution cut-off are 37 and 82 



AU for the iBaraffe et al. 



were 30 and 65 AU). For the 
to 104 and 234 AU. 



(120(331 models, respectiv ely (without the mass correction, these 



Fortney et al 



(120081 ) models these move from 83 and 182 AU 



3.6. Ilda LinI (120041 ) Core Accretion Formation Models 



As in 



Nielsen et al. 



evolution models of ilda fc Lin 



(12008 ) , we 



;urn to the giant planet formation and dynamical 



(120041 ). which predict the final state of giant planets, mass 
and semi-major axis, following the core accretion scenario. We extract 200-300 planets from 
their Fig. 12, and use these masses and semi-major axes in our Monte Carlo simulations. 
We plot the predicted number of planets detected from these models in Fig. [161 with target 
stars divided by binarity. Even wit h our 118 target stars, without removing close binaries. 



not accounting for stellar mass (the 



Ida fc Lin 



star), and using the planet luminosities of the 



(2004) mode 



Baraffe et al 



s were run with a 1 Mq primary 



plan et for each. In 



Nielsen et al. 



mm m odels, the Monte Carlo 
simulations show that for each of t he three cases of Ilda fc Li J ( 2004 ). we would expect to 



(120081 ). we could only exclude the cases 



de tect about 

of Ilda k LinI (|2004| ) A, B, and C at 45%, 49%, and 50% confidence, respectively, with 



the expanded t arget star sam 



63% (using the 



j le her e these rejection levels only b ecome 45%, 59%, an d 



Burrows et al 



When using the models of 



2003) models to be consistent with 



Baraffe et al 



are 38%, 58%, and 62%. The models of 



Nielsen et al. 



mm ). 



(120031). th e limit s for the sample of this paper 



Ida fc LinI (|2004J ) predict very few giant planets 
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in long-period orbits: fewer than 20% of the giant planets predicted by t 
beyon d 10 AU, while our 68% confidence limit on the upper limit of the 



lese models are 



Gumming et al. 



( 120081 ) power-law, 23 AU, gives 30% of giant plane ts in orbits 



current limits can neither confirm nor rule out the 



jeyon d 10 AU. However, our 



Ida fc LinI (12004 ) populations. 



4. Discussion 



Overall, the conclusions from this work are largely similar to those of iNielsen et al. 



(120081 ). that extrasolar giant planets are rare at large separations around Su n-like or less 



mass ive stars. Even with the most pessimistic models for planet NIR fluxes (IFortney et al. 



20081 ). and weighting against the M stars (which provide the most favorable contrasts for 
finding planets), we find that at 95% confidence, fewer than 20% of solar-mass stars can 
have a planet more massive than 4 Mj^p in an orbit between 123 and 218 AU. Also, a 



p ower-law model 
of ICumming et al. 



'o r the semi-major axis distribution of giant planets following the results 

104 AU at 68% confi dence, and 234 AU at 



(120081 ) must have a cut-off o : 



70 confidence, again using the 
mass correction. 



Fortney et al. 



(120081 ) models and the 



Johnson et al. 



mm 



It is worth noting t hat there are ad d itiona l models of giant planet distributions beyond 



what we consider here. 



Gumming et al. 



(120081 ) found a good fit to distributions of close-in 



giant planets from radial velocity work using a single power-law for the distribution of 
semi-major axis. It is also possible, however, that while this is a good fit for giant planets 
within ~5 AU, it may not hold for planets at larger separations; perhaps a broken power 
law or some other distribution governs the population of giant planets in long period orbits. 
If all giant planets are formed beyond the snow line, perhaps the final distribution differs 
for planets that migrate inwards and those that remain beyond the snow line. Alternatively, 
there may be multiple methods of planet formation, such as the core accretion scenario 
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and the disk instability model (e.g. Ilda fc LinI (120041 ) vs. iBosd (120071 )). and the frequency 
with which each occurs is a function of distance from the star. Another possibility is that 
mass and semi-major axis are not independent distributions, which should become testable 
as the number of known exoplanets increases. Additionally, moving across spectral type 



may not only change the frequency of giant planets, but a. 



and semi-major axis. The analysis of iCumming et al. 



so their distributions of mass 



(120081 ) relied exclusively on solar-type 



(FGK) planet hosts, with the number of planets orbiting M stars being too small to draw 
any conclusions about a possible dependence of planet distributions on spectral type. These 
issues cannot be well addressed with further null results; they require a large number of 
detected planets at intermediate (~5-20 AU) and large (>20 AU) separations to make 
statistically significant statements on long-period giant planet populations. 

Since this paper was first submitted, the discovery of several planet candidates, via 
direct imaging, was ann ounced; planets we re detected ar ound the three st ars, all of A 



spectral type, I 



(iLagrange et al 



Johnson et al. 



R 879 9 (IMarois et al 



20081 ) . Fomalhaut (IKalas et al 



20081), and (3 Pic 



20091) . These exciting discoveries are consistent with the predictions of 



(120071 ). as even though planets within 100 AU (like the three found around 



HR 8799) should be easier to detect around lower-mass stars (where the self-luminosity 
of the star is smaller), similar planets have not been found around stars of solar mass or 
smaller. However, each one of these planetary systems is different from the other. Moreover, 
all are around more massive stars than were analyzed here. Since the full survey papers 
for each of these discoveries have not yet been published, it is difficult to incorporate these 
results into our analysis of Sun-like (and less massive) stars. 
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Conclusions 



We have used Monte Carlo simulations to examine the null result from three direct 



imaging surveys (IMasciadri et al 



2005 



Biller et al 



2OO7I : 



Lafreniere et al. 



20071 ) to set 



constraints on the population of extrasolar giant planets. We use three commo nly cited sets 



of planet models (IBurrows et al. 



2003 



Baraffe et al. 



2003 



Fortney et al. 



2008|) in order to 



reach conclusions as broad as possible. Doubling the sample size, as expected, increased 
the strength of our null results. However, including bette r modeling for g iant p lanets-using 



the stellar mass dependence of giant plan et frequency o f 



core-accretion based lu minosity mode l s of 



constraints reported in 



Nielsen et al. 



Fortney et al. 



Johnson et al. 



fl2007h . and the 



(l2008l )-have actually loosened the 



(I2OO8I ). There is still some uncertainty, however, in 



which if any of these models of planet luminosity is correct; likely the truth may fall in 
between that of the optimistic "hot start" models and the somewhat pessimistic "core 
accretion" model. 



With the COND models of iBaraffe et al.l (120031 ) . a planet more massive than 4 Mj^p 



is found around 20% or less of FGKM stars in or bits between 8. 1 and 911 AU, at 68 



confidence. These limits become 7.4 to 863 AU forjBurrows et al. 



for the models of 



Fortney et al. 



fl2003h . and 25 to 557 AU 



(I2OO8I ). At 95% confidence, 4 Mjup (and larger) planets are 



found around fewer than 20% of stars b e tween 22 and 507 AU, 21 an d 479 AU, and 82 an d 



276 AU for the models of 



Baraffe et al 



(120031 ) ■ 



Burrows et al. 



(120031 ) ■ and 



Fortney et al. 



(l2008l ). respectively. 



Using the power law distribution of 



Cumming et al 



(l2008l ). with index -0.61, the upper 



cut-off for the distribution of giant planets is foun d at 30 and 6 5 AU, at the 68% and 95% 



confidence leve l s, resp ectively, using the models of 



Baraffe et al. 



Burrows et al. 



(120031) . With the models o f 



(120031 ) these limits become 28 and 56 AU, and with the 



Fortney et al. 



models they are 83 and 182 AU. 
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When we apply the iJohnson et aLl (j2007f ) dependency of planet fraction on stellar 
mass, the M stars in our sample (where we achieve the greatest sensitivity to planets), 
are weighted down to account for their decreased likelihood of hosting a giant planet. As 
a result, the i mprov ed n ull results cited above r etreat to levels similar to those cited in 



Nielsen et al. 



( 120081 ) and 



Lafreniere et al 



( 120071 ) ■ Given our results, fewer than 20% of 



solar-type sta r s hav e a >4 Mj^p planet between 13 and 849 AU at 68% confidence with t he 



Baraffe et al 



(l2003h mode: 



and 41 and 504 AU for the 
of 



Baraffe et al. 



mm 



s (also 13 and 805 A U for the models of 



Fortnev et a 



Burrows et al. 



Burrows et al 



mm, 



(I20081) m odels). At 95% confid ence, for the models 



(120031), and 



Fortney et al. 



(120081), fewer than 20% of 



1 Mq stars have 4 Mjup planets between 30 and 466 AU, 30 and 440 AU, and 123 and 218 
AU, respectively. 



Applying the 



Johnson et al. 



(120071 ) results to the 



Gumming et al. 



(l2008h model for 



semi-maior a^is dis tribution, giant planets cannot exist beyond 37 and 82 AU for the 



Baraffe et al 



(120031 ) models at 68% a nd 95% confidence. The 68% and 95% confidence 



figures become 3 6 and 82 AU for the .Baraffe et al. 



for the models of 



Fortney et al. 



mm models and 104 and 234 AU 



Johnson et al. 



(120071 ) dependence 



(120081 ). In general, the 
of planet fraction on stellar mass makes direct imaging planet searches more difficult, as 
the stars most likely to harbor giant planets are also the most luminous, giving extreme 
contrast ratios between star and planet that impede planet detection. 

We note that while the constraints on giant planet populations from this and other 
work have, for the first time, reached the equivalent of extrasolar "Kuiper Belts," there is 
still a gap (~5 - ~30 AU) between these results for FGKM stars and those of radial velocity 
surveys, which focus more on the inner solar system. Delving into this unprobed region from 
the direct imaging side can be achieved two ways: increasing sensitivity to planets at small 
separations (achievable with dedicated planet finders using "extreme" adaptive optics, e.g. 
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GPI ( IGraham et al. 



20081 )). or with large-scale 



20071 ) and VLT-SPHERE fiBoccaletti et al.l 
surveys to increase the sample size of target stars, such as the 500 hour Near-Infrared 
Coronagraphic Imager (NICI) survey at Gemini South (Chun et al. 2008). Our technique 
can be applied to the results from any direct imaging survey for giant exoplanets, requiring 
only the target list and achieved contrast curves. By building up the statistics of null 
results, it will be possible to more directly focus direct imaging efforts on where planets 
are most likely to exist, and create a fuller picture of the distribution of extrasolar giant 
planets. Additionally, such an analysis helps to put the survey into context with respect 
to previous work. There is also no limitation based on observation wavelength, even when 
target stars are observed by multiple surveys, since simulated planets are advanced in 
their orbits and compared to each contrast curve. As such, it will be interesting in future 
work to consider th e r esults from L and M ban d surveys currently b eing conducted (e.g. 



Kasper et al. 



fl2007f )l 



Baraffe et al. 



Burrows et al 



(120031 ) predict planets with 



(120031) and 

significantly lower contrasts to their parent stars at these longer wavelengths, and so the 
inclusion of results from such surveys could strengthen our null results, especially at large 



separations. 
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paper, collected 
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Tom the direc t imaging planet 
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Lafreniere et al. 



(120071 ) (Triangles). Tabled] gives other properties of these stars, and Table [2] provides details 
on how the individual ages were determined. The median target star is a 160 Myr Kl star 
at 24 pc. The size and color of the plotting symbols corresponds to the spectral of each 
target star. The top legend gives the conversion between size and color of the plotting 
symbol and spectral type: the color scheme follows the visible spectrum, with early-type 



stars represented by large dark purple symbols, while late-type stars are small red symbols. 



-43- 



-i — I — I — I I I I 



-1 — I — I — I I I I 



U 
Cl 

o 



o 

CD 

o 



10 E 



Boroffe 

15.7 Mj„pO O 

10.5 Mj„pC ^ 

7.53 M,„„o ^ 



20 - 



30 



40 




10 



100 1000 
Age (Myr) 



Fig. 2. 

by the 



A plot of t he ag e and H magnitude o 



Baraffe et al. 



(120031) and 



Burrows et al 



plan ets, for different masses, as predicted 



(I2OO3I ) models, represented by the thin blue 



lines and the thick red lines, respectively. The diamonds and circles are the H magnitudes 
given by the models themselves, while the lines show the interpolation and extrapolations 
beyond these poin ts that we use when assigning H magnitudes to the simulated planets. The 



COND models of 



Baraffe et al. 



(120031 ) required very little extrapolation to f ill the range of 



param eter space shown here, while far more extrapolation is required for the 



Burrows et al. 



(I2OO3I ) models, especially at young ages and small masses. 
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'ig. 3 . — As with Fig. [21 a plot of the predicted fluxes of extrasolar planets of 



Baraffe et al. 



( 20031 ). again represente d by blue lines and 



core accretion models of 



Fortney et al 



open diamonds, this time plotted against the 



(I2OO8I ). the red lines with filled circles. In order to fill 



the 
the 



jarameter spa c e of p lanet mass and stellar age we consider, it is necessary to extrapolate 



Fortney et al. 



(120081 ) H magnitudes beyond the grid points of the models themselves, 



especially at larger ages and smaller masses. 
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Fig. 4. — The upper limit on planet fraction {fp, the fraction of stars with a planet of a given 
mass and semi-major axis, see Equation [1]), at the 95% (bl ue, thin lines) and 68 % (red, thick 



Baraffe et al 



(120031). With 95% 



lines) confidence levels, using the theoretical models of 
confidence, we can say that less than 1 in 20 stars has a planet more massive than 8 Mj^p 
between 50 and 160 AU (constrained by the solid blue, thin curve). We plot a horizontal 
fiducial bar (again, with a thick red line and thin blue line) at 4 Mjup, intersecting the 
fp <20% contour at both 68% confidence (outer contour, thick red line) and 95% (inner 
contour, thin blue line). Hence, the horizontal line at the bottom right of the figure suggests 
no more than 1 in 5 stars would have a planet more massive than 4 Mjup from 8.1 to 911 
AU at the 68% confidence level, and between 22 and 507 AU at the 95% confidence level. 
Known radial velocity planets are shown as filled circles for comparison. 
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Fig. 5. 



As with Fig, ffl the u pper hmit on planet fraction only now using the theoretical 



models of Burrows et al 



(120031 ). The overall shape of the graph is quite similar, so with 95% 
confidence, we can place an upper limit on planet fraction of 5% for planets larger than 8 
Mjup with semi-major axis between 55 and 130 AU. As before, radial velocity planets are 
plotted as solid circles. The fiducial fp <20% limits for 4 Mj^p are between 7.4 and 863 AU 
at 68% confidence, and 21 to 479 for the 95% confidence level. 
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Fig. 6. — The same as Fig. [Hand Fig. [5l this time using the models of 



to giv e the upper hmit on planet fraction. Overall, the theoretical models of iFortney et al. 



Fortnev et al. 


(2008) 


dels of iFortnev et al. 



(120081 ) are more pessimistic as to NIR fluxes o f planets when compared to the hot-start 



models (IBurrows et al. 



2003 



Baraffe et al. 



20031). With these models, from 82 to 276 AU, 



less than 20% of stars can have a planet above 4 Mj^p, at the 95% confidence level, and 
between 25 and 557 AU at 68% confidence. 
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Fig. 7. — Twelve models for the semi-maio r axis distribution of extrasolar planets, using the 



planet luminosity models of 



Baraffe et al 



(120031 ) . with power law indices of a = -1, -0.61, 
and -0.25, and upper cut-offs (the limit up to which there are planets, but beyond which 
planets no longer appear) of 10, 20, 40, and 80 AU. The solid purple line gives the histogram 
of known radial velocity planets, the horizontal and diagonal green lines give different values 
of the power law index, and the red vertical lines mark the upper cut-offs. The vertical 
black dashed line at 2.5 AU gives the approximate upper limit to which the radial velocity 
survey is complete to planets. The percentages at each intersection of power law and upper 



cut-off show the confidence with which that model 



' dN 
■■ da 



oc a" for a < a, 



cut-off, 



and 



dN 
da 







for a > acut-off) can be rejected. For example, a planet population with dN/da ~ a and 



an outer cutoff o 



(ICumming et al 



10 A U is ruled out at 5.3% confidence. For the power law of index -0.61 



20081 ). at 95% confidence the upper cut-off must be less than 65 AU, which 



would fall between the dashed and dot-dashed vertical lines of this graph. 
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Fig. 8. — Contours showing the confidence with which we can exclude models of the semi- 
major axis distribution of extrasolar giant planets of the form ^ oc a", with an upper cut-off 



beyond which there are no longer pl anets, using the r node 



power law index of -0.61 as given by 



Gumming et al. 



s of 



Baraffe et al. 



fl2003f ). The 



(120081 ) is marked with a dotted line. 



The jags in the contours are due to binaries being removed as we move up in power law 
cut-off (binary target stars are pulled once the considered semi-major axis cut-off reaches 
one-fifth the binary separation). The pronounced jag between 50 and 55 AU corresponds to 
the binary M-dwarfs TWA 8A and TWA 8B (21 pc, 10 Myr) being removed from the sample, 
indicating the strong effect a few M stars have on our results. For the power law of index 
-0.61, the 68% and 9 5% confidence l evels for rej ection of this model are at 30 and 65 AU. 



Similar plots for the 



Burrows et al 



fl2003h and 



Fortney et al 



(120081 ) models are available 



in our supplement, available at http://exoplanet.as.arizona.edu/~lclose/exoplanet2.html 



which also contains versions of all future plots for these additional two sets of models. The 



68% and 95% confid ence levels are 28 and 56 AU (Burrows et al. 



( jFortney et al. 



20031) and 83 and 182 AU 



20081 ). 
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Fig. 9. — Our linear fit to the dependence of tlie likeliliood a target star lias of hosting 
a close-in, giant extrasolar planet as a function of stellar m ass. The histogram shown is 



the metallicity - correc ted histogram of 



Johnson et al 



Johnson et al. 



(120071 ) (their Fig. 6). As noted by 



(120071 ). the probability for the high-mass bin is likely underestimated, so 
future work may show an even greater boost for the value of high-mass target stars. For the 
target stars considered in this work, 35 are in the low-mass bin, 78 are in the medium-mass 



bin, and 5 are in the bin for the highest masses. 
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Fig. 10. — Contours giving the upper limit on planet fraction around a. 



Baraffe et al 



1 the stars of M 



(120031). Comparing to 



spectral type in the three surveys, using the models of 
Fig. m which considered stars of all spectral types, the behavior of the contours at small 
semi-major axis is roughly the same, while the outer edge and depth of the upper limit are 
limited by the reduced sample size (only 18 of the 118 target stars are M stars). For 68% 
confidence, fewer than 1 in 5 stars have a planet more massive than 4 Mj„p between 9.0 and 



207 AU. For the models of 



AU for the iFortney et al. 



Burrows et al. 



(12008 



(12003h this range is 8.3 to 213 AU, and is 43 to 88 



) models. 
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Fig. 11. — The upper limit on planet fraction, using only the FGK stars in o ur survey (as 



well a s the single A star in the survey, HD 172555 A), using the models of iBaraffe et al. 



(120031 ). The shapes of the contours and the behavior at large semi- major axes are roughly 
the same as in Fig. HJ when all stars were considered, but without the M stars and their 
favorable contrasts at small separations, the small-period planets are much less accessible. 
The 20% contours, at the 68% and 95% confidence levels, for planets more massive than 4 
M j„p, are foun d , betw een 25 and 856 AU and between 38 and 469 AU, respectively. For the 



Burrows et al. 



(I2OO3I ) model s, the 68% an d 95% limit ranges are between 25 and 807 AU, 



(I2OO8I ) models, the 95% confidence 20% contour 



and 40 and 440 AU; for the iFortney et al. 
never reaches 4 Mjup, but the 68% confidence range is from 59 to 497 AU. 
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Fig. 12. — The upper limit on planet fraction for stars of 1 Mq (/p,i.o), with c onstraints from 



target stars of higher and lower mass stars weighted according to a fit to the iJohnson et al. 
(20o3) de pendence on st e llar m ass of the frequency of radial velocity planets, using the 



models of 



Burrows et al 



(I2OO3I ). The plot is similar to that of Fig. |H which weighted all 
stars equally, but the contours shrink slightly (mainly on the low separation side of the plot) 
as the M stars are now effectively given less weight. The 20% confidence level for planets 
more massive than 4 Mj„p are between 13 and 8 49 AU at 68% confidence, and between 30 



Burrows et al. 



and 466 AU for the 95% confidence level. For the 
are 13 to 805 AU, and 30 to 440 AU, while for the models of 
are between 41 and 504 AU, and 123 and 218 AU. 



fl2003h models, t 



Fortney et al 



lese ranges 



20081) the limits 



Baraffe, Mass Correction to 0.5 M© 




Fig. 13. — The upper limit on planet fraction , this time normaliz ing; to stars of a half solar 



mass (fp,o.5); or about MO, using the models of iBaraffe et al.l (120031) . The constraints become 



stronger, as expected, as our assumption going into this calculation is that lower mass stars 
are less likely to have planets overall. With this set of assumptions, the dearth of giant, large- 
separation planets around M stars is made quite clear. The possibility of lower mass, inner 
planets around M stars (and indeed, planets like those in our own solar system) remains, 
however. Fewer than 20% of M stars can have planets more massive than 4 Mj^p between 



9.0 and 1070 AU at 68% confidence, and 23-605 AU at the 95% confidence 



limits are 8.3 to 1016 AU and from 22 to 573 AU for the iBur 



well as 26 to 656 AU and 71 to 341 AU for the 



Fortney et al. 



•ows 



et al 



fl2008h 



evel. These 



(120031) models, as 
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Fig. 14. — As with Fig.[7l power-law models for the semi-major axis distribution of extrasolar 
planets, and the confidence wit h which we c a n rule out these models, using the results of 



our survey an d 



Johnson et al 



the models of 



Baraffe et al. 



(120031 ). This time, we utilize the results of 



20071 ) to appropriately give additional weight to high mass stars, which are 



more likely to harbor giant planets. 
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Fig. 15. — Contours showing the confidence with which we can exclude a model for the 
distribution of semi-major axis of extrasolar giant plane ts given by ^ oc a" up to some upper 



cut-off. This figur e shows the results fo r the models of 



Baraffe et al. 



(120031), with the stellar 



(120071 ) to account for the dependence of like 



mass correction of Johnson et al. 
giant planets upon the mass of the parent star. The upper cut-off for the 



ihood of finding 



Gumming et al. 



(120081 ) power law of index -0.61 (as marked by the dotted 



Burrows et al 



i ne) is 37 AU at the 68% confidence 



(120031 ) mode 



level, and 82 at 95% confidence. For the 
confidence limits are at 36 and 82 AU, and at 104 and 234 AU for the 
models. 



s, the 68% and 95% 



Fortney et al 



(120081) 
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Fig. 16. — The number of planets we'd expect to detect, as a function of the number of stars 
in our survey. Target stars are divided into bins: one for single stars, and binaries divided by 
separation; within each bin the best targets are pla ced to the left of the graph, so they're "ob- 



Baraffe et al 



served first" in this manner. Using the models of 

for stellar mass effects or removing binaries, the core accretion models of 
predic t detecting between 0.5 and 1 planets w i th the combination of the 



(120031 ) . a nd not accounting 



Ida fc Linl (120041) 



(l2005h 



Biller et al. 



(l2007h . and 



Lafreniere et al. 



20071 ) surveys. The 



Masciadri et al. 



Ida fc Lin 



(|200J) mod- 



els A, B, and C (green curves) can only be excluded with 38%, 58%, and 62% confidence, 
respectively. 
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Table 1. Target Stars 



Target 




RA^ 




Dcc^ 


Distance (pc)^ 


Sp. Type 


Age (Myr) 




H^ 


Ks^ 


Obs. Mode^ 


Biller et al. (2007) 


HIP 1481 


00 


18 


26, 


,1 


-63 28 39.0 


40.95 


F8/G0V 


30 


7.46 


6, 


.25 


6.15 


VLT SDI 


HD 8558 


01 


23 


21, 


,2 


-57 28 50.7 


49.29 


G6V 


30 


8.54 


6, 


.95 


6.85 


VLT SDI 


HD 9054 


01 


28 


08, 


.7 


-52 38 19.2 


37.15 


KIV 


30 


9.35 


6, 


.94 


6.83 


VLT SDI 


HIP 9141 


01 


57 


48, 


.9 


-21 54 05.0 


42.35 


G3/G5V 


30 


8.11 


6, 


.55 


6.47 


VLT SDI 


BD+05 378 


02 


41 


25, 


.9 


-|-05 59 18.4 


40.54 


MO 


12 


10.20 


7, 


.23 


7.07 


VLT SDI 


HD 17925 


02 


52 


32, 


. 1 


-12 46 11.0 


10.38 


KIV 


200 


6.05 


4, 


.23 


4.17 


VLT SDI/GDPS 


Eps Eri 


03 


32 


55, 


.8 


-09 27 29.7 


3.22 


K2V 


1100 


3.73 


1, 


.88 


1.78 


VLT SDI/GDPS 


V577 Per A 


03 


33 


13 


.5 


-f46 15 26.5 


33.77 


G5IV/V 


70 


8.35 


6, 


.46 


6.37 


MMT SDI 


GJ 174 


04 


41 


18, 


.9 


+20 54 05.4 


13.49 


K3V 


280 


7.98 


5, 


.31 


5.15 


VLT SDI 


GJ 182 


04 


59 


34, 


.8 


+01 47 00.7 


26.67 


MlVc 


12 


10.10 


6, 


.45 


6.26 


VLT SDI/Ks/GDPS 


HIP 23309 


05 


00 


47, 


,1 


-57 15 25.5 


26.26 


MO/1 


12 


10.09 


6, 


.43 


6.24 


VLT SDI/Ks 


AB Dor 


05 


28 


44, 


.8 


-65 26 54.9 


14.94 


K2Vk 


70 


6.93 


4, 


.84 


4.69 


VLT SDI 


UY Pic 


05 


36 


56, 


.8 


-47 57 52.9 


23.87 


KOV 


70 


7.95 


5, 


.93 


5.81 


VLT SDI 


AO Men 


06 


18 


28, 


,2 


-72 02 41.4 


38.48 


K6/7 


12 


10.99 


6, 


.98 


6.81 


VLT SDI/Ks 


HIP 30030 


06 


19 


08, 


.1 


-03 26 20.0 


52.36 


GOV 


30 


8.00 


6, 


.59 


6.55 


MMT SDI 


HIP 30034 


06 


19 


12, 


.9 


-58 03 16.0 


45.52 


K2V 


30 


9.10 


7, 


.09 


6.98 


VLT SDI 


HD 45270 


06 


22 


30, 


.9 


-60 13 07.1 


23.50 


GIV 


70 


6.50 


5, 


.16 


5.05 


VLT SDI 


HD 48189 A 


06 


38 


00, 


.4 


-61 32 00.2 


21.67 


G1/G2V 


70 


6.15 


4, 


.75 


4.54 


VLT SDI 


piOl UMa 


08 


39 


11, 


.7 


-f 65 01 15.3 


14.27 


G1.5V 


200 


5.63 


4, 


.28 


4.17 


MMT SDI/GDPS 


HD 81040 


09 


23 


47, 


.1 


-|-20 21 52.0 


32.56 


GOV 


2500 


7.74 


6, 


.27 


6.16 


MMT SDI 


LQ Hya 


09 


32 


25, 


.6 


-11 11 04.7 


18.34 


KOV 


13 


7.82 


5, 


.60 


5.45 


MMT/VLT SDI/Ks/GDPS 


DX Leo 


09 


32 


43, 


.7 


-|-26 59 18.7 


17.75 


KOV 


200 


7.01 


5, 


.24 


5.12 


MMT/VLT SDI/GDPS 


HD 92945 


10 


43 


28, 


.3 


-29 03 51.4 


21.57 


KIV 


70 


7.76 


5, 


.77 


5.66 


VLT SDI/GDPS 


GJ 417 


11 


12 


32, 


.4 


-|-35 48 50.7 


21.72 


GOV 


200 


6.41 


5, 


.02 


4.96 


MMT SDI/GDPS 


TWA 14 


11 


13 


26, 


.5 


-45 23 43.0 


46.00^ 


MO 


10 


13.00 


8, 


.73 


8.49 


VLT SDI 


TWA 25 


12 


15 


30, 


.8 


-39 48 42.0 


44.00^ 


MO 


10 


1 1.40 


7, 


.50 


7.31 


VLT SDI 


RXJ1224.8-7503 


12 


24 


47, 


.3 


-75 03 09.4 


24.17 


K2 


16 


10.51 


7, 


.84 


7.71 


VLT SDI 


HD 114613 


13 


12 


03, 


.2 


-37 48 10.9 


20.48 


G3V 


8800 


4.85 


3, 


.35 


3.30 


VLT SDI 


HD 128311 


14 


36 


00, 


,6 


+09 44 47.5 


16.57 


KO 


630 


7.51 


5, 


.30 


5.14 


MMT SDI 


EK Dra 


14 


39 


00, 


,2 


+64 17 30.0 


33.94 


GO 


70 


7.60 


6, 


.01 


5.91 


MMT SDI/GDPS 


HD 135363 


15 


07 


56, 


,3 


+76 12 02.7 


29.44 


G5V 


3 


8.72 


6, 


.33 


6.19 


MMT SDI/GDPS 


KW Lup 


15 


45 


47, 


,6 


-30 20 55.7 


40.92 


K2V 


2 


9.37 


6, 


.64 


6.46 


VLT SDI 


HD 155555 AB 


17 


17 


25, 


,5 


-66 57 04.0 


30.03 


G5IV 


12 


7.20 


4, 


.91 


4.70 


VLT SDI/Ks 


HD 155555 


17 


17 


27, 


,7 


-66 57 00.0 


30.03 


M4.5 


12 


12.70 


7, 


.92 


7.63 


VLT SDI/Ks 


HD 166435 


18 


09 


21, 


,4 


+29 57 06.2 


25.24 


GO 


110 


6.85 


5, 


.39 


5.32 


MMT SDI 


HD 172555 A^ 


18 


45 


26, 


,9 


-64 52 16.5 


29.23 


A5IV/V 


12 


4.80 


4, 


.25 


4.30 


VLT SDI 


CD -64 1208 


18 


45 


37, 


,0 


-64 51 44.6 


34.21 


K7 


12 


10.12 


6, 


.32 


6.10 


VLT SDI/Ks 


HD 181321 


19 


21 


29, 


,8 


-34 59 00.5 


20.86 


G1/G2V 


160 


6.48 


5, 


.05 


4.93 


VLT SDI 


HD 186704 


19 


45 


57, 


,3 


+04 14 54.6 


30.26 


GO 


210 


7.03 


5, 


.62 


5.52 


MMT SDI 


GJ 799B 


20 


41 


51, 


.1 


-32 26 09.0 


10.22 


M4.5c 


12 


11.00 


5, 


.20 


-99.00 


VLT SDI/Ks 


GJ 799A 


20 


41 


51, 


,2 


-32 26 06.6 


10.22 


M4.5c 


12 


10.25 


5, 


.20 


4.94 


VLT SDI/Ks 


GJ 803 


20 


45 


09, 


.5 


-31 20 27.1 


9.94 


MOVc 


12 


8.81 


4, 


.83 


4.53 


VLT SDI/Ks/GDPS 


HD 201091 


21 


06 


53, 


,9 


+38 44 57.9 


3.48 


K5Ve 


2000 


5.21 


2, 


.54 


2.25 


MMT SDI 


Eps Indi A 


22 


03 


21, 


.7 


-56 47 09.5 


3.63 


K5Ve 


4000 


4.69 


2, 


.35 


2.24 


VLT SDI 


GJ 862 


22 


29 


15, 


,2 


-30 01 06.4 


15.45 


K5V 


6300 


7.65 


5, 


.28 


5.11 


VLT SDI 
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Table 1 — Continued 



Target 




RAl 


Dec^ 


Distance (pe)^ 


Sp. Type 


Age (Myr) 


vi 


h3 


Ks3 


Obs. Mode* 


HIP 112312 A 




22 44 57.8 


-33 15 01.0 


23.61 


M4e 


12 


12.20 


7.15 


6.93 


VLT SDI 


HD 224228 




23 56 10.7 


-39 03 08.4 


22.08 


K3V 


70 


8.20 


6.01 


5.91 


VLT SDI 


Masciadri et al. 


(2005) 




















HIP 2729 




00 34 51.2 


-61 54 58 


45.91 


K5V 


30 


9.56 


6.72 


6.53 


VLT Ks 


BD +2 1729 




07 39 23.0 


02 11 01 


14.87 


K7 


200 


9.82 


6.09 


5.87 


VLT H/GDPS 


TWA 6 




10 18 28.8 


-31 50 02 


77.00^ 


K7 


10 


11.62 


8.18 


8.04 


VLT Ks 


BD +1 2447 




10 28 55.5 


00 50 28 


7.23 


M2 


70 


9.63 


5.61 


5.31 


VLT H/GDPS 


TWA 8A 




11 32 41.5 


-26 51 55 


21.00^ 


M2 


10 


12.10 


7.66 


7.43 


VLT Ks 


TWA 8B 




11 32 41.5 


-26 51 55 


21.00^ 


MS 


10 


15.20 


9.28 


9.01 


VLT Ks 


TWA 9A 




11 48 24.2 


-37 28 49 


50.33 


K5 


10 


11.26 


8.03 


7.85 


VLT Ks 


TWA 9B 




11 48 24.2 


-37 28 49 


50.33 


Ml 


10 


14.10 


9.38 


9.15 


VLT Ks 


SAO 252852 




14 42 28.1 


-64 58 43 


16.40^ 


K5V 


200 


8.47 


5.69 


5.51 


VLT H 


V343 Nor 




15 38 57.6 


-57 42 27 


39.76 


KOV 


12 


8.14 


5.99 


5.85 


VLT Ks 


PZ Tel 




18 53 05.9 


-50 10 50 


49.65 


KOVp 


12 


8.42 


6.49 


6.37 


VLT Ks 


BD-17 6128 




20 56 02.7 


-17 10 54 


47.70 


K7 


12 


10.60 


7.25 


7.04 


VLT Ks 


Lafreniere et al. 


(2007) 




















HD 166 




00 06 36.7839 


-1-29 01 17.406 


13.70 


KOV 


200 


6.13 


4.63 


4.31 


GDPS 


HD 691 




00 11 22.4380 


-1-30 26 58.470 


34.10 


KOV 


260 


7.96 


6.26 


6.18 


GDPS 


HD 1405 




00 18 20.890 


+30 57 22.23 


30.60 


K2V 


70 


8.60 


6.51 


6.39 


GDPS 


HD 5996 




01 02 57.2224 


-H69 13 37.415 


25.80 


G5V 


440 


7.67 


5.98 


5.90 


GDPS 


HD 9540 




01 33 15.8087 


-24 10 40.662 


19.50 


KOV 


2900 


6.96 


5.27 


5.16 


GDPS 


HD 10008 




01 37 35.4661 


-06 45 37.525 


23.60 


G5V 


200 


7.66 


5.90 


5.75 


GDPS 


HD 14802 




02 22 32.5468 


-23 48 58.774 


21.90 


GOV 


5200 


5.19 


3.71 


3.74 


GDPS 


HD 16765 




02 41 13.9985 


-00 41 44.351 


21.60 


F7IV 


290 


5.71 


4.64 


4.51 


GDPS 


HD 17190 




02 46 15.2071 


-1-25 38 59.636 


25.70 


KlIV 


4300 


7.81 


6.00 


5.87 


GDPS 


HD 17382 




02 48 09.1429 


-1-27 04 07.075 


22.40 


KIV 


430 


7.62 


5.69 


5.61 


GDPS 


HD 18803 




03 02 26.0271 


-1-26 36 33.263 


21.20 


G8V 


4400 


6.72 


5.02 


4.95 


GDPS 


HD 19994 




03 12 46.4365 


-01 11 45.964 


22.40 


F8V 


6200 


5.06 


3.77 


3.75 


GDPS 


HD 20367 




03 17 40.0461 


-1-31 07 37.372 


27.10 


GOV 


380 


6.41 


5.12 


5.04 


GDPS 


2E 759 




03 20 49.50 


-19 16 10.0 


27.00 


K7V 


200 


10.26 


7.66 


7.53 


GDPS 


HIP 17695 




03 47 23.3451 


-01 58 19.927 


16.30 


M3e 


70 


11.59 


7.17 


6.93 


GDPS 


HD 25457 




04 02 36.7449 


-00 16 08.123 


19.20 


F5V 


70 


5.38 


4.34 


4.18 


GDPS 


HD 283750 




04 36 48.2425 


-1-27 07 55.897 


17.90 


K2 


300 


8.42 


5.40 


5.24 


GDPS 


HD 30652 




04 49 50.4106 


-1-06 57 40.592 


8.00 


F6V 


4500 


3.19 


1.76 


1.60 


GDPS 


HD 75332 




08 50 32.2234 


-1-33 17 06.189 


28.70 


F7V 


270 


6.22 


5.03 


4.96 


GDPS 


HD 77407 




09 03 27.0820 


-1-37 50 27.520 


30.10 


GO 


120 


7.10 


5.53 


5.44 


GDPS 


HD 78141 




09 07 18.0765 


-1-22 52 21.566 


21.40 


KO 


270 


7.99 


5.92 


5.78 


GDPS 


HD 90905 




10 29 42.2296 


-1-01 29 28.025 


31.60 


GOV 


230 


6.90 


5.60 


5.52 


GDPS 


HD 91901 




10 36 30.7915 


-13 50 35.817 


31.60 


K2V 


1000 


8.75 


6.64 


6.57 


GDPS 


HD 93528 




10 47 31.1553 


-22 20 52.927 


34.90 


KIV 


310 


8.36 


6.56 


6.51 


GDPS 


HIP 53020 




10 50 52.0645 


-1-06 48 29.336 


5.60 


M4 


200 


11.66 


6.71 


6.37 


GDPS 


HD 96064 




11 04 41.4733 


-04 13 15.924 


24.60 


G8V 


250 


8.41 


5.90 


5.80 


GDPS 


HD 102392 




11 47 03.8343 


-11 49 26.573 


24.60 


K4.5V 


3400 


9.05 


6.36 


6.19 


GDPS 


HD 105631 




12 09 37.2563 


-1-40 15 07.399 


24.30 


KOV 


1500 


8.26 


5.70 


5.60 


GDPS 


HD 107146 




12 19 06.5015 


-1-16 32 53.869 


28.50 


G2V 


190 


7.07 


5.61 


5.54 


GDPS 


HD 108767 B 




12 29 50.908 


-16 31 14.99 


26.90 


K2V 


140 


8.51 


6.37 


6.24 


GDPS 



Table 1 — Continued 



Target RA^ Doc^ Distance (pc)^ Sp. Typo Age (Myr) Ks'^ Obs. Mode* 



HD 109085 


12 


32 


04 


.2270 


-16 11 45.627 


18. 


,20 


F2V 


100 


4.31 


3, 


,37 


3, 


,37 


GDPS 


BD +60 1417 


12 


43 


33 


.2724 


+60 00 52.656 


17. 


,70 


KO 


270 


9.40 


7, 


,36 


7, 


,29 


GDPS 


HD 111395 


12 


48 


47, 


.0484 


+24 50 24.813 


17. 


,20 


G5V 


1000 


6.31 


4, 


,70 


4, 


,64 


GDPS 


HD 113449 


13 


03 


49, 


.6555 


-05 09 42.524 


22. 


10 


KIV 


70 


7.69 


5, 


,67 


5, 


,51 


GDPS 


HD 116956 


13 


25 


45 


.5321 


+56 58 13.776 


21. 


,90 


G9V 


710 


7.29 


5, 


,48 


5, 


,41 


GDPS 


HD 118100 


13 


34 


43 


.2057 


-08 20 31.333 


19. 


,80 


K4.5V 


280 


9.31 


6, 


,31 


6, 


,12 


GDPS 


HD 124106 


14 


11 


46 


.1709 


-12 36 42.358 


23. 


10 


KIV 


1700 


7.92 


5, 


,95 


5, 


.86 


GDPS 


HD 130004 


14 


45 


24, 


.1821 


+ 13 50 46.734 


19. 


,50 


K2.5V 


5100 


7.60 


5, 


,67 


5, 


,61 


GDPS 


HD 130322 


14 


47 


32, 


.7269 


-00 16 53.314 


29. 


,80 


KOIII 


2900 


8.05 


6, 


,32 


6, 


,23 


GDPS 


HD 130948 


14 


50 


15 


.8112 


+23 54 42.639 


17. 


,90 


G2V 


420 


5.88 


4, 


,69 


4, 


,46 


GDPS 


HD 139813 


15 


29 


23 


.5924 


+80 27 00.961 


21. 


,70 


G5 


270 


7.31 


5, 


,56 


5, 


,45 


GDPS 


HD 141272 


15 


48 


09, 


.4630 


+01 34 18.262 


21. 


,30 


G9V 


280 


7.44 


5, 


,61 


5, 


,50 


GDPS 


HIP 81084 


16 


33 


41 


.6081 


-09 33 11.954 


31. 


,93 


K9Vkce 


70 


11.29 


7, 


78 


7, 


,55 


GDPS 


HD 160934 


17 


38 


39, 


.6261 


+61 14 16.125 


24. 


,54 


K7 


70 


10.18 


7, 


,00 


6, 


,81 


GDPS 


HD 166181 


18 08 


16.030 


+29 41 28.12 


32. 


,58 


G5V 


60 


7.70 


5, 


,61 


5, 


,61 


GDPS 


HD 167605 


18 


09 


55, 


.5001 


+69 40 49.788 


30. 


,96 


K2V 


500 


8.60 


6, 


,45 


6, 


,33 


GDPS 


HD 187748 


19 


48 


15, 


.4478 


+59 25 22.446 


28. 


,37 


GO 


140 


6.66 


5, 


,32 


5, 


,26 


GDPS 


HD 201651 


21 


06 


56 


.3893 


+69 40 28.548 


32. 


,84 


KO 


6800 


8.20 


6, 


,41 


6, 


,34 


GDPS 


HD 202575 


21 


16 


32 


.4674 


+09 23 37.772 


16. 


17 


K3V 


700 


7.91 


5, 


,53 


5, 


,39 


GDPS 


HIP 106231 


21 


31 


01 


.7137 


+23 20 07.374 


25. 


,06 


K3Vkc 


70 


9.24 


6, 


,52 


6, 


,38 


GDPS 


HD 206860 


21 


44 


31, 


.3299 


+ 14 46 18.981 


18. 


,39 


GOVCH-0.5 


200 


6.00 


4, 


,60 


4, 


,56 


GDPS 


HD 208313 


21 


54 


45 


.0401 


+32 19 42.851 


20. 


,32 


K2V 


6400 


7.78 


5, 


,68 


5, 


.59 


GDPS 


V383 Lac 


22 


20 


07 


.0258 


+49 30 11.763 


10. 


,68 


KO 


40 


8.57 


6, 


,58 


6, 


,51 


GDPS 


HD 213845 


22 


34 


41 


.6369 


-20 42 29.577 


22. 


,74 


F5V 


200 


5.20 


4, 


,27 


4, 


,33 


GDPS 


HIP 114066 


23 


06 


04, 


.8428 


+63 55 34.359 


24. 


,94 


MO 


70 


10.87 


7, 


17 


6, 


,98 


GDPS 


HD 220140 


23 


19 


26 


.6320 


+79 00 12.666 


19. 


74 


K2Vk 


85 


7.73 


5, 


,51 


5, 


,40 


GDPS 


HD 221503 


23 


32 


49 


.3999 


-16 50 44.307 


13. 


,95 


K6Vk 


550 


8.60 


5, 


,61 


5, 


,47 


GDPS 


HIP 117410 


23 


48 


25 


.6931 


-12 59 14.849 


27. 


,06 


K5Vkc 


55 


9.57 


6, 


,49 


6, 


,29 


GDPS 



"'"from the CDS Simbad service 

^derived from the Hipparcos survey Iperryman et alj IiQQtI 



^from the 2MASS Survcv lcutri et alJ hoO^ 

'^In cases where target stars were observed by multiple surveys, the star is listed only in the first section of this table where it appears, 
cither in the lBiller et alJ j^OoJ) or llVIasciadri et alJ l200^ section, with Observing Mode given as "VLT SDI/Ks" or "VLT H/GDPS," for 
example. 

^Distance from lsong et alJ J2OO3I ) 

'^As this is the only star in our sample earlier than F2, we consider this work to be a survey of FGKM stars. 
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Table 2. Age Determination for Target Stars 



Target Sp. Type* LiEW(mas)* Li Age (Myr) ^'HK* R'/fit Age Group Membership"^ Group Age^ Adopted Age + + + 



Biller et al. (2007) 


HIP 1481 


FS/GOV^ 


129'^ 


100 


-4.360^ 


221 


Tuc/Hor 


30 


30 


HD 8558 


G6V^ 


205^ 


13 






Tuc/Hor 


30 


30 


HD 9054 


KIV^ 


170^ 


160 


-4.236 


100 


Tuc/Hor 


30 


30 


HIP 9141 


G3/G5V''' 


181^ 


13 






Tuc/Hor 


30 


30 


BD+05 378 


MO^ 










/3 Pic 


12 


12 


HD 17925 


KIV^ 


194 


50 


-4.357 


216 


Her/Lyr 


200 


200 


Eps Eri 


K2V^'^ 






-4.598^ 


1129"'' 






1100 


V577 Per A 


G5IV/V"'""'" 


219 


3 






AB Dor 


70 


70 


GJ 174 


K3V^^ 


45^ 


280 


-4.066 








280 


GJ 182 


MlVe"'"'* 


280"'"'^ 


12 










12 


HIP 23309 


MO/1"'"'^ 


294-'-'^ 


12 


-3.893 




/3 Pic 


12 


12 


AB Dor 




267^ 


10 


-3.880^ 


<50 


AB Dor 


70 


70 


UY Pic 


KOV"'-® 


263^ 


10 


-4.234'^ 


78 


AB Dor 


70 


70 


AO Men 


K6/7-'-^ 


357-'-^ 


6 


-3.755*^ 




(3 Pic 


12 


12 


HIP 30030 


GOV"-'-^ 


219^ 


2 






Tuc/Hor 


30 


30 


HIP 30034 


K2V^ 










Tuc/Hor 


30 


30 


HD 45270 


GlV^ 


149'^ 


90 


-4.378 


254 


AB Dor 


70 


70 


HD 48189 A 


G1/G2V^ 


145^ 


25 


-4.268^ 


105 


AB Dor 


70 


70 


piOl UMa 


Gl.SV^*^ 


135^ 


100 


-4.400^"'^ 


300 






200 


HD 81040 


GOV^*^ 


2^22 


2500 










2500 


LQ Hya 


KOV^*^ 


247^ 


13 










13 


DX Leo 


KOV^*^ 


180^ 


100 


-4.234 


78 


Her/Lyr 


200 


200 


HD 92945 


KIV^*^ 


138^ 


160 


-4.393^ 


285 


AB Dor 


70 


70 


GJ 417 


GOV^"^ 


76 


250 


-4.368^"^ 


235 


Her/Lyr 


200 


200 


TWA 14 


MO^^ 


600^^ 


8 






TW Hya 


10 


10 


TWA 25 


MO^ 


494 


10 






TW Hya 


10 


10 


RXJ1224.8-7503 


K2^^ 


250^"^ 


16 










16 


HD 114613 


G3V^^ 


100 


400 


-5.118 


7900 






8800 


HD 128311 


KO^" 






-4.489^^ 


565 






630 


EK Dra 


GO^" 


212^ 


2 


-4.106^^ 


<50 


AB Dor 


70 


70 


HD 135363 




220^ 


3 










3 


KW Lup 


K2V^^ 


430^^ 


2 










2 


HD 155555 AB 


G5IV1^ 


205^ 


6 


-3.965^ 


<50 


Pic 


12 


12 


HD 155555 C 


M4.5^^ 










Pic 


12 


12 


HD 166435 


GO^^ 






-4.270^^ 


107 






110 


HD 172555 A 


A5IV/V^ 










Pic 


12 


12 


CD -64 1208 




580^^ 


5 






Pic 


12 


12 


HD 181321 


G1/G2V^^ 


131^ 


79 


-4.372^ 


243 






160 


HD 186704 


GO^^ 






-4.350^^ 


205 






210 


GJ 799B 


M4.5e^^ 










Pic 


12 


12 


GJ 799A 


M4.5e^'^ 










Pic 


12 


12 


GJ 803 




51« 


30 






Pic 


12 


12 


HD 201091 


K5Ve34 






-4.704^^ 


2029+ 






2000 


Eps Indi A 


K5Ve3^ 






-4.851^ 


3964+ 






4000 


GJ 862 


K5V3^ 






-4.983^ 


6280+ 






6300 
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Table 2 — Continued 



Target 




Sp. Typo* 


Li EW (mas)* 


Li Age (Myr) 




R'/jj^ Age"'"^ Group Membership^ 


Group Age^ 


Adopted Age + + 


HIP 112312 A 




M4e^ 








13 Pic 


12 


12 


HD 224228 




K3V^* 


53* 


630 


-4.468'' 


AB Dor 


70 


70 


Masciadri et al. 


(2005) 
















HIP 2729 




K5V2 








Tue/Hor 


30 


30 


BD +2 1729 












Her/Lyr 


200 


200 


TWA 6 






560^'^ 


3 




TW Hya 


10 


10 


BD +1 2447 




M2^^ 








AB Dor 


70 


70 


TWA 8A 




M2^^ 


530^° 


3 




TW Hya 


10 


10 


TWA 8B 




M6^^ 


560^'^ 


3 




TW Hya 


10 


10 


TWA 9A 




K5^^ 


460^'' 


3 




TW Hya 


10 


10 


TWA 9B 




Ml^^ 


480^^ 


3 




TW Hya 


10 


10 


SAO 252852 




K5V^'' 








Her/Lyr 


200 


200 


V343 Nor 




KOV^ 


300^^ 


5 


-4.159'' 


40 (3 Pic 


12 


12 


PZ Tel 




KOVpl^ 


267^* 


20 


-3.780* 


<50 13 Pic 


12 


12 


BD-17 6128 






400*" 


3 




13 Pic 


12 


12 


Lafreniere et al. 


(2007) 
















HD 166 




K0V*1 


7442 


290 


-4.458^^ 


460 Her/Lyr 


200 


200 


HD 691 




KOV*"^ 


110* 


260 


-4.380^^ 


260 




260 


HD 1405 




K2V*'* 


271*'^ 






AB Dor 


70 


70 


HD 5996 




G5V*'' 






-4.454^^ 


440 




440 


HD 9540 




KOV'' 






-4.774'^ 


2900 




2900 


HD 10008 




G5V''^ 


103** 


280 


-4.530^^ 


740 Her/Lyr 


200 


200 


HD 14802 




GOVl^ 


51*^ 


4000 


-4.985'' 


6300 




5200 


HD 16765 




F7IV''° 


73^* 


270 


-4.400^^ 


300 




290 


HD 17190 




KlIV^^ 






-4.870^1 


4300 




4300 


HD 17382 




KIV^^ 






-4.450^^ 


430 




430 


HD 18803 




G8V^2 






-4.880^1 


4400 




4400 


HD 19994 




F8V^^ 


12^4 


8000 


-4.880^^ 


4400 




6200 


HD 20367 




GOV^^ 


113* 


150 


-4.500^1 


610 




380 


2E 759 






63^'' 


260 




Her/Lyr 


200 


200 


HIP 17695 




M3e^^ 








AB Dor 


70 


70 


HD 25457 




F5V^'' 


gj60 


80 


-4.390^^ 


280 AB Dor 


70 


70 


HD 283750 




K28I 


33* 


300 


-4.057^^ 






300 


HD 30652 






15" 


7500 


-4.650^^ 


1500 




4500 


HD 75332 






125* 


50 


-4.470^1 


500 




270 


HD 77407 




GO^^ 


162*= 


50 


-4.340^^ 


190 




120 


HD 78141 




K08* 


107* 


270 








270 


HD 90905 




GOV^^ 


136* 


80 


-4.430^^ 


370 




230 


HD 91901 




K2V'' 


766 


1000 








1000 


HD 93528 




KIV^'' 


100* 


260 


-4.424'' 


360 




310 


HIP 53020 




M4S7 








Her/Lyr 


200 


200 


HD 96064 




GSV**** 


114* 


250 


-4.373^^ 


250 




250 


HD 102392 




K4.5Vl" 






-4.811'' 


3400+ 




3400 


HD 105631 










-4.650^^ 


1500 




1500 


HD 107146 




G2V™ 


125* 


180 


-4.340^1 


190 




190 


HD 108767 B 




K2V''^ 


175^2 


140 








140 



-I- 
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Table 2 — Continued 



Target 


Sp. Type* 


Li EW (mas)* 


Li Age (Myr) 


R'hk' 


H'ffJf Age++ 


Group Membership* 


Group Age* 


Adopted Age + + 


HD 109085 


F2Vl^ 


3^73 


100 










100 


BD +60 1417 




96** 


270 










270 


HD 111395 








-4.580""- 


1000 






1000 


HD 113449 


KIV^'' 


142* 


200 


-4.340*^ 


190 


AB Dor 


70 


70 


HD 116956 




Si"-* 


1000 


_4.447l3 


420 






710 


HD 118100 


K4.5V''8 


25*^ 


280 


-4.090*^ 








280 


HD 124106 


KIV^'' 






-4.675" 


1700 






1700 


HD 130004 


K2.5V''8 






-4.919*3 


5100+ 






5100 


HD 130322 


KOIII^^ 






-4.780"* 


2900 






2900 


HD 130948 


G2V52 


116* 


230 


-4.500"* 


610 






420 


HD 139813 


G5™ 


119* 


230 


-4.400"* 


300 






270 


HD 141272 








-4.390"* 


280 






280 


HIP 81084 








-4.210'^ 




AB Dor 


70 


70 


HD 160934 




40^9 


280 






AB Dor 


70 


70 


HD 166181 


GSV**" 


186* 


60 










60 


HD 167605 


K2V81 


14" 


500 










500 


HD 187748 




114* 


140 










140 


HD 201651 


KO*'' 






-5.010"* 


6800 






6800 


HD 202575 


K3V''8 






-4.522*3 


700+ 






700 


HIP 106231 


K3Vkc'=8 


140* 


180 


-3.906*3 




AB Dor 


70 


70 


HD 206860 


GOVCH-0.5^^ 


110*" 


190 


-4.400'^ 


300 


Her /Lyr 


200 


200 


HD 208313 


K2V'5** 






-4.987*3 


6400+ 






6400 


V383 Lac 




259* 


40 










40 


HD 213845 








-4.547'^ 


830 


Her /Lyr 


200 


200 


HIP 114066 


MO'^8 










AB Dor 


70 


70 


HD 220140 


K2Vk''8 


218* 


85 


-4.074*3 








85 


HD 221503 


K6Vk^'' 






-4.486" 


550+ 






550 


HIP 117410 


K5Vkcl'' 






-4.194" 


55+ 






55 



*Group Membership for TWA, /3 P ic, Tuc/Hor, and AB Dor fro m IZuckerman fc Sonj 12004^ . Her/L yr from 
Ages from lzuckerman fc Sonj hoO'h (TWA, Pic, and Tuc/Horl . iNielsen et al 



r— T. 



* Measurement References: 2 : [HoiU^^^^owle^ _^^7; 



Houk fc S mith-M oord lll98 



i200j ), 12: 



Leaton i 



Gray et al.| i2006bn. 18: 



: Page] 



Bidelmad 



A 1 1951I 1. 



Wichmann et al.l 



13: 



2009 1. 23: 

I982I). 29: iRandich et al.l 



iGliese fc JahreisJiiggill. 35: 
jiggsl). 40: iMathioudakis e~al. 



J1993I 



), 3: 



Gray et al.l i2003a 



197? ), 19: 



Gaidos 



I), 30: 



et al 



Waite et al 



1 20051) 



20051 ) (AB Dor), and [L6pez-Sant 
1994- 5: 



Lopez-Santiago et al.| j200l 



e^^antiagf 
iago et al.| 



2003) . 9: IZuckerman 8z Son. 
) 



Outispoto et al 



Gliese fc Jahreifl 



1199' 



36: 



H_^^^et_al, 

J 2004I) . U 
I995I1 



Zuckerman et al.l i2001b 



3 11979I- 



I), 31: 



Torres et al.H2000l1. 6: 



et al 



1196 



1995), 20: Montes et al.li 2001bl). 2 1: 



Favata et al.l h997l). 16: 



JTll: 



1 20061) (Her/Lyr) 

12006 J). 



Gray et al 



Christian Mathioudakiq 

1 



Zuckerman et al 



, .Wright et al 

Song et abl J2003l ). 27: 



45: 
50: 



Montes et al 



IQQcfl. 56: 



3 



Cowlevl Jl97£), 51 



koOlaj). 46: 
Ineardl 



Fleming et al 



[oswalt et al.l il988l) . 62: 

lelma nl jigssl). 



Rufener fc BarthoMil I1982I). 



Neuhauser fc Brandnei 
1943 )^_37^ 



Helmer et al.l h983l'l. 47 



J Il96ll] . 



32: 



Egge 



1 I2OO4I). 



22: 



Alcala et al.l I 199; ) , 28: 



|198£ ), 57: Favata 
ay et al.l 2003bll 



Mallik et al.l 




^^^^a^^^J I 1995 ) , 58 : ^^^nenzelle^^^^^ 
rrj J1969I ' 



19561). 52: iHarlan fc Taylor! il970bl)y 53: iHerbig fc Spald 



all jl997l' 



Zboril et al 



h, 48: 



1991 

38: ISoderblom et al 
19971). 43: 



19961 



L^gez^^^ia^o_e^^L|_^00^ 



33: 



Sozzetti et al 



19981). 39: 



119851 



Nesterov et al.l 



Il99i 



in j I I955I' 



59: 



I ). 64: ISchwope et al.l i2000l ). 65: 



19731 



70: 



75: lUpgren 



I), 44: lAmbruster et al 
\) 



I). 54: [Israel ian et al 



Harlan &z Taylod il97o'a|v" 

Staro nl ]l97d ). 76: |pye et al 



I 197; ), 60 , 



Mora et al 



I1995I) . 



197 



1 



PasQuini et al 



20041). 55 




Sato fc Kuiil 



Lambert fc Reddvl i2004l 



Strassmeier et al 



2001 ). 72: Pallavicini et al 



77: I Fan et al 



78: 



2000 



1995) 
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J 19951 ). 79: Izuckerman et alj J2OO4I '}, 80: lEggeiJ jlQej) . 81: Istocke et al] ]l99l|) . 82: I Chen et al] kooj) . 83: leowyer etal\ 1x99^ ) 

"'"In general, wo have only determined Ca R'/f ages for stars with spectral types Kl or earlier, but in the case of these K2-K6 stars, we have 
only the R'/f/f nieasurement on which to rely for age determination. The calibration of Mt. Wilson S-index to h K stars {B-V 1.1 mag) 

has not been wcll-dcfincd jNoyes et al.| jl984|) ; spccifieally the photosphcric subtraction), and hence applying a R'/fj^- vs. age relation for K5 stars 



is unlikely to yield accurate ages. 

"''"'"Using Eq. 3 of llVIamaiek &z Hillenbrandl ]2008l ) to convert h K hito age 

"'""'""'"In general, ages derived from lithium and/or calcium alone are likely accurate to within a factor of '^2 
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Table 3. Summary of Results. 



Target Stars Mass Correction* Confidence Level ^araff^^^^^ (^00^) ^urrow^^^^^^ ^00^) ^Wtne^^e^^^ ^00^ 



Completeness plots: semi-major axis range with fp < 20% for M > 4 Mj. 



All 


None 




68% 


8.1 - 911 AU 




7.4 - 863 AU 


25 - 557 AU 


All 


None 




95% 


22 - 507 AU 




21 - 479 AU 


82 - 276 AU 


M stars 


None 




68% 


9.0 - 207 AU 




8.3 - 213 AU 


43 - 88 AU 


M stars 


None 




95% 










FGK stars 


None 




68% 


25 - 856 AU 




25 - 807 AU 


59 - 497 AU 


FGK stars 


None 




95% 


38 - 469 AU 




40 - 440 AU 




All 


1 M0 




68% 


13 - 849 AU 




13 - 805 AU 


41 - 504 AU 


All 


1 Mq 




95% 


30 - 466 AU 




30 - 440 AU 


123 - 218 AU 


All 


0.5 Mq 




68% 


9.0 - 1070 AU 




8.3 - 1016 AU 


26 - 656 AU 


All 


0.5 Mq 




95% 


23 - 605 AU 




22 - 573 AU 


71 - 341 AU 




Upper cut-off on power law 


distribution for semi-major 


axis with index -0.61 




All 


None 




68%, 


30 AU 




28 AU 


83 AU 


All 


None 




95% 


65 AU 




56 AU 


182 AU 


All 


Yes 




68% 


37 AU 




36 AU 


104 AU 


All 


Yos 




95% 


82 AU 




82 AU 


234 AU 


The "Mass Correction" column 


refers 


to whether or 


not theljohnson et al.l j 


2001 


\) result, that more mass: 


ive stars are more likely 


to harbor giant planets, is used to 


wcigfit 


the target stars by stellar mass. For 


the 


completeness plots, this 


eorreetion is either not 


applied (None) or 


set to a specific 


stellar 


mass, to dete: 


rmine the upper limit on the frequency of giant planets around stars of that 


mass. For the lim. 


its on tire upper 


cut-off 


on power law 


■ distributions, the eorreetion is either applied (Yes) 


or not (None). 
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Table 4. Binaries 



Target 


Sep {■■)! 


Sep. (AU)1 


Reference 


Companion Type 


Biller ct al. f2007) 


HIP 9141 


0.15 


6.38 


Biller et al. (2007) 


mid-G 


\ \j 1 1 i cl Jr\. 


7 


230 


Pounds et al. fl993) 


MO 


AB Dor 


9 (Ba/Bb) 


134 (Ba/Bb) 


Close et al. f2005) 


Binary M stars 


AB Dor 


0.15 (C) 


2.24 (C) 


Close et al. f 20051 


Very low-mass M Star 


HIP 30034 


5.5 


250 


Chauvin et al. (2005) 


Planet/Brown Dwarf 


HD 48189 A 


0.76 (B) 


16.5 


Fabricius & Makarov f20001 


K star 


HD 48189 A 


0.14 


3.03 


Biller et al. (2007) 


K star 


DX Leo 


65 


1200 


Lowrance et al. ('20051 


M5.5 


EK Dra 


SB 


SB 


Metchev & Hillenbrand (2004) 


M2 


HD 135363 


0.26 


7.65 


Biller et al. f20071 


late K/early M 


HD 155555 AB 


SB (AB) 


SB (AB) 


Bennett et al. (1967) 


G5 and KO SB 


HD 155555 AB 


18 (C) 


1060 (C) 


Zuckerman et al. r200la1 


Target Star 155555 C, M4.5 


HD 172555 A 


71 


2100 


Simon & Drake (1993) 


Target Star CD -64 1208, K7 


HD 186704 


13 


380 


Aitken & Doolittle (19321 


early M 


GJ 799 A 


3.6 


36 


Wilson (19541 


Target Star GJ 799B, M4.5 


HD 201091 


16 


55 


Baize (19501 


K5 


Eps Indi A 


400 


1500 


McCauKhrean et al. (2004) 


Binary Brown Dwarf 


HIP 112312 


100 


2400 


Sone et al. (20021 


M4.5 


Masciadri et al. ('2005) 


TWA 8A 


13 


270 


Javawardhana et al. (1999) 


Target Star TWA 8B, M5 


TWA 9A 


9 


576 


Javawardhana et al. (19991 


Target Star TWA 9B, Ml 


SAO 252852 


15.7 


260 


Poveda et al. (1994) 


HD 128898, Ap 


V343 Nor 


10 


432 


Sons et al. (2003) 


M4.5 


BD-17 6128 


2 


100 


Neuhauser et al. (2002) 


M2 


Lafreniere et al. (2007) 


HD 14802 


0.47 


10 


Lafreniere et al. (2007) 


K6 


HD 16765 


4.14 


89 


Holden (1977) 


~K 


HD 17382 


20.3 


456 


Leoine & Shara (2005) 


M4.5 


HD 19994 


~5 


~100 


Hale (1994) 


M3V 


HD283750 


124 


2220 


HolberE et al. (2002) 


White Dwarf 


HD 77407 


1.7 


50 


Musrauer et al. (20041 


~M 


HD 93528 


234 


8200 


Perrvman et al. (19971 


HIP 52776, K4.52 


HD 96064 


11.47 


283 AU 


LioDincott & MacDowall (19791 


NLTT 26194, M3 


HD 102392 


1.13 


28 


Lafreniere et al. (20071 


~M 


HD 108767 B 


23.7 


639 


Gould & Chaname (20041 


AOIV 


HD 130948 


2.64 


47 


Potter et al. (20021 


Binary brown dwarfs 


HD 139813 


31.5 


683 


Stephenson (19601 


GO 


HD 141272 


17.8 


350 


Eisenbeiss et al. (20071 


MS 


HD 160934 


SB 


SB 


Hormuth et al. (20071 


early M, a — 4.5 AU 


HD 160934 


8.69 


213 


Lowrance et al. (2005) 


~M 


HD 166181 


SB 


SB 


Nadal et al. (1974) 


1.8 day orbit 


HD 166181 


0.102 


3 


Lafreniere et al. (20071 


~K 


HD 167605 


1.2 


37 


Arribas et al. (19981 


M4V 


HD 206860 


43.2 


795 


Luhman et al. (2007b1 


T dwarf 


HD 213845 


6.09 


139 


Lafreniere et al. (20071 


late M 
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Table 4 — Continued 



Target 


Sep C')! 


Sep. (AU)1 


Reference 


Companion Type 


HD 220140 


10.9 


214 


Lowranee et al. (2005) 


mid M 


HD 220140 


963 


19000 


Makarov et al. (2007) 


~M 


HD 221503 


339 


4700 


Gould & Chaname (2004) 


binary M stars 


HIP 117410 


1.84 


50 


Rossiter (1955) 


early M 



"'"SB indicates a spectroscopic binary 

^Thcse stars have Hipparcos proper motion and parallax within errors, and similar values 
of calcium R' hk {-4.424 and -4.451 for HD 93528 and HIP 52776, respectively). 



